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Abstract 

The main goal in this paper is to point out that quantity ||Vi?|p(p) 
on a harmonic space can not be determined by the spectra of local 
geodesic spheres or balls, therefore the main results of |AM-S| (quoted 
in the title) are wrong. My strong interest in this theorem is mo- 
tivated by the fact that it contradicts some of my isospectrality ex- 
amples constructed on geodesic spheres and balls of certain harmonic 
manifolds. The authors overlooked that the Lichnerowicz identity is 
not determined by the given spectral data, and so is the final crucial 
equation obtained by eliminating with the Lichnerowicz identity. In 
short, the above theorem has falsely been established by spectrally un- 
determined identities which can not be computed (determined) by the 
spectra of local geodesic spheres. More complicated spectrally unde- 
termined functions cause the problems in case of local geodesic balls. 
I describe also a strong physical argument which clearly explains why 
the manifolds appearing in my examples are isospectral. 

It must be pointed out, however, that a very remarkable new idea, 
namely, the asymptotic expansion of the heat invariants ak {p, r) de- 
fined on geodesic spheres, Sp{r), is introduced in the paper. It can be 
used for developing both geometric uncertainty theory and global vs. 
local spectral investigations. Among my contributions to this develop- 
ing field is the following statement: Average volumes, J vol {Bp (r)) dp 
resp. / vol{Sp{r))dp, of geodesic balls resp. spheres are generically 
not determined by the spectra of compact Riemann manifolds. This 
theorem interestingly contrasts the statement asserting that the vol- 
ume of the whole manifold can be determined in terms of this global 
spectrum. In other words: Music written for local drums can not be 
played, in general, on the whole manifold. 
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1 A short history of harmonic manifolds. 



Several equivalent definitions of harmonic manifolds are known in the lit- 
erature. By one of them, they are the Riemannian manifolds yielding the 
local mean value theorem of harmonic functions. By an other characteri- 
zation, they are the manifolds where the density function ^ det{gij){q) = 
Op{q)/r^~^ , defined on local normal coordinate systems about an arbitrary 
point p G M", is radial, that is, it depends just on the distance r{p, q). One 
can proof then that density function 6p{r) does not depend on the reference 
point p either. 

The latter characterization clearly shows that all locally two-point ho- 
mogeneous manifolds are harmonic. The converse direction concerned, Lich- 
nerowicz (1946) showed that harmonic manifolds satisfying dimM^ = n < 4 
are exactly the locally 2-point homogeneous manifolds. The question if this 
statement is true also in higher dimensions, became known as Lichnerowicz 
conjecture concerning harmonic manifolds. In 1990, this problem was affir- 
matively answered on compact manifolds having finite fundamental groups, 
by this author in [Szl] . Due to the fact that all harmonic manifolds are Ein- 
stein, they have constant scalar curvature, which is obviously strictly positive 
in the above compact cases. By Myers' and Cheeger-Gromoll's theorems, 
this affirmative answer extends also onto complete harmonic manifolds sat- 
isfying Seal > |Sz2) . However, extension onto non-compact manifolds 
satisfying Seal < defied every effort used by this author. 

Then, in 1992, without knowing anything about the Lichnerowicz con- 
jecture, Damek and Ricci |DR] accidentally discovered infinitely many non- 
compact and locally non-symmetric harmonic manifolds. Their statement 
says that the natural invariant metrics on the solvable extensions, SH^j^'^\ 

of Heisenberg-type groups, H^'^'^\ are harmonic. Such manifolds exist for 
any I > 0. Later we investigate families defined by fixed values (a -|- b) and 
I. If I ^ 3 mod 4, the metrics in such a family are isometric. However, 
if / 7^ 3 mod 4, two metrics in the same family are locally non-isometric, 
unless one of them is associated with (a, 6) and the other with (6, a). Par- 
ticularly interesting are the families defined for / = 3. In fact, the metric 
on SH^°'~^^'^^ ~ SH^'""^^^ is symmetric while the others are locally non- 
symmetric. 

Due to this discovery, the investigations were intensified both on com- 
pact and non-compact harmonic manifolds. Since the Damek-Ricci exam- 
ples do not have compact factors, the question arises if the conjecture can 
be established on all compact manifolds and not just on those satisfying 
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Seal > 0. By using entropy theory of geodesic flows, this question has been 
positively answered on closed manifolds having negative sectional curvature, 
by G. Besson, G. Courtois and S. Gallot |BCG] . in 1995. Since manifolds 
of strictly negative constant scalar curvature can have sectional curvature 
of changing sign, this question is still open on closed manifolds belonging 
to this category. In his very recent paper |Kn] . Knieper confirms the Lich- 
nerowicz conjecture for all compact harmonic manifolds without focal points 
or with Gromov hyperbolic fundamental groups. These results strongly sug- 
gests that the conjecture must be true on all manifolds that have compact 
Riemannian factors. The Damek-Ricci spaces have also been very intensely 
investigated in the past two decades. In this short list of achievements let 
it be mentioned just J. Heber's result |H], who, in 2006, proved that a ho- 
mogeneous harmonic space is either a model space or a Damek-Ricci space. 
The existence of inhomogeneous harmonic manifolds is an other important 
open question in this field. 

By this author, spectral investigations of harmonic manifolds were initi- 
ated. Since then, this field has been developed into several different direc- 
tions. This paper can focus mostly on the investigations of this author. They 
inevitably involve also the Heisenberg type groups whose solvable extensions 
are the Damek-Ricci spaces. Several different isospectrality examples have 
been constructed. The isospectral domains arise always on the members of a 
family H^'^'^^ (resp. SH^"''''^) defined by the same (a + b) and The AM&S- 
paper refers to the examples constructed in [Sz41 ISz5j . where the isospec- 
trality is confirmed for any pair of geodesic balls resp. spheres having the 
same radius r. The AM&S-theorem contradicts the examples constructed 
on the particular family SH^"''''\ where the metric on SH^""^^'^^ ~ SH^'^^^^ 
is symmetric while the others in the isospectrality family are locally non- 
symmetric. What makes this situation more serious is that, in 2005, Hagen 
Fiirstenau [Fj discovered mathematical difficulties in the construction of the 
intertwining operator introduced in |Sz4l ISz5j . But, soon thereafter, this 
author did come up with the solution of this problem |Sz6t ISz8] , which was 
publicly announced also at a conference held at CUNY, in February, 2006 
[Sz7j . Since then, I found also a strong physical argument explaining why 
the isospectrality stated in |Sz41 ISz5] must be true. Namely, it turns out, 
that the Laplacian on the investigated manifolds can be identified with the 
Hamilton operators of elementary particle systems and the isospectrality is 
equivalent to the C-symmetry obeyed by these physical systems. This solu- 
tion (together with explaining the difficulties in |Sz4t ISz5] ) are described in 
the last section of this paper. 
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2 Technicalities on harmonic manifolds. 



Along a unit speed geodesic c(r), the Jacobian endomorphism field ylc(r)(-) 
acting on the Jacobi subspace Vc(^r) -L c(r) = u{r) is defined by the equations 

K{r) + Rc{r)Ac{r) = 0, ^c(O) = 0, -4^(0) = (1) 

where i?c(r)(-) = Ru{-) ■= R{-,u)u is the Jacobian curvature along c(r). 
Its normalized version is denoted by Ac(j.) = ^^^(r). The power series of 
invariants Tr A^^^^ can directly be computed by these formulas. For instance, 
for k = 1 we have: 




+-^«.(o)^c,(o) + 2i?c,(o)4,(o) + ^^^(o)) + 0{r'), 

where Cp(r) indicates that it starts out from p = Cp(0) and R'^ (o)(') ~ 
Kpi-) = (V„pi?)(.,np)up, etc. 

Among these invariants the density Oc[o){c{r)) = det ^^(j,) plays the most 
important role in this paper. Its power expansion is usually computed by 
Legendre's recursion formulas. Before providing the corresponding formulas 
computed by this direct method, we have to point out that some of the basic 
objects introduced here are different from those appearing in j AM-S] or |Be] . 
As opposed to these texts, the Laplacian is a negative semi-definite operator 
and the curvature is defined by R{X, Y)Z = VxVyZ— Vy V[x,y]-2^ in 

all of my papers quoted later. In order to keep the exposition simple, in this 
section we keep the AriasMarco-Schiith-definition, R{x,y)z = —S/x^yZ + 
S/yVxZ + V[x,y]Z, for the curvature, but will return back to our definition 
in the rest part of this paper. There is a difference also regarding 9 in 
this paper which is associated with detA in |AM-S| . It should also be 
mentioned that some of the constants appear in |Be) incorrectly not just in 
those formulas which have been corrected in |AM-Sj but also in some other 
expressions which are quoted in this paper (these errors and their effects on 
my computations were pointed out to me by D. Schiith). In terms of these 
quantities we have: 

Theorem 2.1. (see IBe'j, Chapter 6). If M is harmonic then there exist 
constants C, H, L ^ such that for all p & M and all u € TpM with \u\ = 1 
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we have: 



Tr(Ru) = C; in particular : 
Tr{R^^) = 0, for all derivatives of order /c G N. 

Tr{RuRu) = H] in particular : 
Tr{RuR'u) = and 
Tr(i?„0 = -Tr{R'M- 
Tr{2,2RuRuRu - 9«) = L. 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



The first five conditions imply that harmonic manifolds are Einstein, 
furthermore, the curvature has constant norm ||-R|p = RahcdR"'^'^'^ ■ Due 
to the Einstein property, these manifolds are analytic and the unique radial 
density 0{r) is uniquely determined by the constants appearing in the Taylor 
expansion. Such a density function does not determine the metric. See, for 
instance, that the above families of Damek-Ricci spaces share the very same 
density function 0{r). An other surprising feature of these examples is that in 
families Sh''^'^^ member SH^"'~^^'^^ is symmetric while the others are locally 
non-symmetric. 

In this paper particular attention is paid to condition ([8]). It should be 
emphasized that this condition declares this curvature-expression to be con- 
stant just for this combination but it says nothing about the terms RuRuRu 
and R'uR'u which can very well be non-constant functions even along a fixed 
geodesic. The same density function can manifest itself very differently. For 
instance, on S ' \ constant L is incorporated completely into Ry^R^^R^ 
but on the other members the same L is distributed both into RuRuRu and 
R'^R'^. Even the proportion of the two terms depends on the members of a 
family. 

By integrating these constants on the unit sphere formed by unit vectors 
u in the tangent space Tp{M) one eliminates direction u and obtains the 
curvature conditions depending just on p G M. This computations result 
the following well known formulas: 



where R and R are defined in an orthonormal basis {ei, . . . , e„} of the tan- 
gent space by the formulas: 



R\\^ = f^{(n + 2)H-C'), (9) 



32{nC'^ + -C\\R\\^ + -R - R) 



27||Vi?|p = n(n2 + 6n + 8)L, (10) 




(11) 



Once again, the harmonicity conditions impose just this single equation 
onto a tensor field having a large number of components. Thus, beyond this 
relation, it provides no further information about the main players; R, R, 
and ||V-R||^; in this formula. In other words, these individual functions are 
out of touch of the harmonicity conditions. Due to Lichnerowicz, there is 
an other independent equation: 

- ^A\\R\\^ = 2C\\R\\^ - R-4R+\\VR\\^ (12) 

which holds true not just on harmonic but on all Einstein manifolds. On 
harmonic manifolds, relation = constant implies (see [Be] , Proposition 

6.68): 

112^-32^ = 27||Vi?|p + Qi, R + AR=\\\/Rf + Q2, (13) 

therefore, R=^\\VRf + Q3, R=^\\VR\\^ + Qi, (14) 
24 yb 

where Qi,Q2tQ3, Q4 are constants which can be expressed in terms of C, H, 
and L. (The second formula on line (llSp appears with an incorrect factor 2 
(before ||Vii:|p) in [Be].) 

It should be pointed out that there is a substantial difference between the 
equations (jlOp and (|12p . The first one is defined by integrating the constant 
function ([8]) on the unit sphere while no such natural integral representation 
exist for the Lichnerowicz identity. This identity is originally established 
by operations defined just on M which are nothing to do with curvature 
quantities such as R'^R'^;...; e. t. c, defined in terms of directions u. This 
detachment from functions defined by radial expansions strongly indicates 
that the Lichnerowicz identity may cause serious difficulties in investigating 
spectrally determined identities obtained by the radial expansion of the heat 
invariants ak{p,r). In fact, it will turn out that this identity is not encoded 
into the spectra of local geodesic spheres, by any means, and its application 
for reformulating spectrally determined identities produces spectrally unde- 
termined ones. This is why the identity established in |AM-S) in the final 
step is not spectrally determined. Therefore, it can not prove the desired 
spectral determinacy stated in the paper. There is also an other explanation 
for this problem which does not use Lichnerowicz identity at all. Although it 
is defensible in case of geodesic spheres, it develops into a major argument 
against the AM&S-statement established by the spectra of local geodesic 
balls. The precise details cirG clS follows. 
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3 The AriasMarco-Schiith local heat invariants. 



The first few heat invariants of a compact manifold S are given by: 

ao{S) = vol{S), ai{S) = - [ scaldvols, (15) 

6 Js 

02(5) = -^ [ {bscal'^ -2\\Ric\\'^ + 2\\R\\'^)dvols. (16) 
360 Js 

They are defined as coefficients in the asymptotic expansion 

00 00 
Tr{exp{tA)) = exp{Xit) ~ (4^t)-'ii"^(5)/2 ^ ak{S)t'' (17) 

1=0 k=0 

where t J, and Aj — )• — 00 denote the eigenvalues of Laplacian A with 
multiplicities. 

Arias-Marco and Schiith furnished the spectral data of local geodesic 
spheres by the asymptotic expansions of radial functions ak{Sp{r)), where 
Sp{r) denotes the geodesic sphere of radius r about p. The first few terms are 
explicitly computed for a2{Sp{r)). The following theorem is a combination of 
Propositions 3.2 and 4.2 of |AM-Sj . however, it is not a literal transcription 
of those statements. In order to prepare the arguments discussed below, this 
transcription focuses also on the eliminating steps (i24]) - (f27|) . These are just 
briefly explained in |AM-Sj without introducing any formulas. 

Theorem 3.1. JAM-S f Consider an n- dimensional harmonic space M with 
constants C,H;L, introduced above. Let p € M and let u be a unit vector 
inTp{M). Then 

ll^^c'^''^''''llLp(r«) =a-4r"^ + a_2r"^ + ao + a2r^ + 0(r^), (18) 
and I \R^^^^^ I = /3_4r-4 + f3_2r~^ + /3o + (32r^ + 0{r^) (19) 

for r 4 0, where the coefficients Ui and f3i for i = —4, —2, 2 are constants 
depending only on n,C and H. Moreover, 

«2(^)=«2 + 3^rr(«) (20) 

4 " 

and /32(^i) = /32 + g J^Tr(i?„oi?(ei,.)i?«ei), (21) 

i=l 

where the coefficients 02 and $2 are constants depending only on n;C; H; L 
and {ei, . . . , e^} is an orthonormal basis of Tp{M). 
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By integrating on the unit sphere S'i(Op) C Tp{M), whose volume is 
denoted by ujn-i, we get: 



[ TriR'Mdu = ^"IwV ,J |Vfi|P(p), (22) 
JsiiOp) n(n + 2)(n + 4) 

„ n 

/ V TriRu o R{ei, .)Ruei)du = (nC^ + 2R{p) - -R{p))^±^. (23) 
Js,(0.)i-1 4 n{n + 2) 



vl i=l 

By using for eliminating terms R{p) and R{p)) from i23\) we get: 
2R{p) - -^R{p) = ^\\VR\\\p) + Q{n,C,H,L), (24) 

where Q is a constant depending only on constants appearing in its argument, 
moreover: 

a2.= -^j a2{u)du = a2 + -—^——-A\VR\\\p), (25) 
^n-i Jsi(Op) 16n(n + 2)(n + 4) 

If .... . 9 ,2, 



P2-= / /32(n)dn = /32 + — — — -|lVi?||^(p), (26) 

^n-i y5i(0p) 32n(n + 2) 

where 0.2 and j32 are constants depending only on n,C,H, and L. Thus the 
scrutinized coefficient is rehorn-ed in the new form: 

The computation proceeds by inserting this asymptotic expansion of 
a2{Sp{r)) into the integral formula /g^,-^-, a2{Sp{r)) ^^°II^Sp(r)) ' where the in- 
tegration is taken with respect to the normalized density '^^0° (g'^ff)"* • It is 
particularly important at this step that the density is constant regarding r, 
that is, it appears as the normalized density on the Euclidean unit sphere. 
Therefore, in the asymptotic expansion of the whole integral formula this 
density has no effect on the expansion. Also note that, due to the uniquely 
determined density function 9{r) on a fixed harmonic manifold (or for a 
whole family SH^""'^^), the constants like C, H, L,...e. t. c. are spectrally 
determined values. After term by term integration of the corresponding co- 
efficient of the asymptotic expansion, the process is finished by eliminating 
R{p) and R{p) by formulas (|14p . Then the final conclusion is this: Since 
the coefficients of the asymptotic expansion as well as the constant values 
are spectrally determined, therefore also ||V-R|p(p) must be spectrally de- 
termined. This statement is considered also regarding the spectra of local 
geodesic balls. This case is discussed in the next section. 
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4 Examining the AM&S-proof. 



4.1 Case of geodesic spheres. 

It must be emphasized again that the concerns brought up against this proof 
appear in the final ehminating step. Ah preceding, highly non-trivial com- 
putations are performed with impressive carefulness. The key problem has 
already been indicated at the introduction of the Lichnerowicz identity. Ac- 
cording to those comments, this equation does not have any natural integral 
representation and its application cuts the computations completely off from 
the spectral data. It seems to be that it is not encoded into the spectra of 
local geodesic spheres at all, thus, by adding it to a local heat invariant, it 
may change it into a quantity which is not local heat invariant any more. 
Next, we show that exactly this is the case. We will prove that this identity 
is really a spectrally undetermined object and when it is added, in the form 
of the above described linear combination with (jlOp . to the local heat invari- 
ant computed in the above theorem, then this elimination process changes 
this spectrally determined invariant to a spectrally undetermined curvature 
expression which can not be produced as linear combination of "true" local 
heat invariants obtained by direct radial expansions of several heat invari- 
ants ak{p,r). The conclusion with such a "fake" local heat invariant gives, 
of course, a false proof for the desired spectral deter minacy of ||V-R||. 

In order to establish these statements, one should define the spectral 
data at a fixed point p G M. It is defined by the linear space of curvature 
expressions obtained by radial expansions of heat invariants ak{p,r). They 
can be computed by integrating the coefficients defined by the radial ex- 
pansions of the heat invariant functions regarding dup. The evaluations of 
these integrals provide the curvature expressions which can be written up 
as identities such that the right sides of these identities are thought to be 
the integral values which are equal to the curvature expressions obtained by 
the AriasMarco-Schiith expansions. On the left sides, the different curva- 
ture symbols can be considered as independent objects (unknowns). These 
curvature terms are defined on M, that is they depend on p and not on the 
directions Up. This identity interpretation of these curvature expressions 
better indicates where these "true" local heat invariants come from. Thus 
the given spectral data at p is a linear space spanned by the identities ob- 
tained by the AriasMarco-Schiith radial expansions of the complete set of 
heat invariants. 

In order to define the complete linear space of local heat invariants at 
p, one should introduce a naturally defined pre-Hilbert norm on this linear 
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space. It is defined by integrating the products Fi{u)F2{u) of infinitesimal 
lieat invariant functions by means of du. By the topological closure associ- 
ated with this pre-Hilbert norm, one can define the spectral-identity- space 
SIS{p), which Hilbert space contains all spectrally determined identities at- 
tached to the given spectral data. The identities belonging to SIS{p) are 
the only ones which are determined by the given spectral data. Also note 
that this spectral data determines just the identities obtained by integrating 
the infinitesimal heat invariant functions, but, it does not determines the 
pre-functions F(u) from which these identities are derived. 

Because of this new interpretation, the identities introduced so far must 
be rewritten in new appropriate forms. This reformulation mostly concerns 
the constants appearing on the left sides of these identities. Also these 
terms must appear as combinations of curvature terms (unknowns) obtained 
by integrating the pre-functions. The values standing on the right sides 
are just formally symbolized by the integral / F{u)dun, where dun denotes 
normalized measure, whose actual evaluation provides the curvature terms, 
depending just on p, on the left sides of these equations. To this end, 
consider the constants C{u),H{u), and L{u) defined in formulas Q-dH]) as 
functions of u which are actually defined by the left sides of those equations. 
Then L on the right side of (fTO|) must be substituted by / L{u)dun while 
and C||-R|p on the left side should be replaced by the curvature expressions 
obtained by actual evaluations of integrals / C^{u)dun and ^ / C{u){{n + 
2)H{u) — C'^{u))dun- These integrals provide linear combinations of terms 
such as Scal\j{p) and ScalM{p)\\R\\^{p)- On non-Einstein manifolds these 
expressions are more complicated involving also the Ricci tensor and some 
other curvature terms. Similar reformulation should be implemented in 
(fT2]l . (f20]l . ([2T]) . (I23])-(126]) and also radial function scat^ in ([THD should be 
expanded. 

It is important to understand that this reformulation is not just a kind 
of fussiness. For instance, identities obtained by evaluating integrals like 
/ C^{u)dun and ^ / C{u){{n + 2)H{u) - C'^{u))dun are not encoded into 
the spectra of local geodesic spheres, meaning that they can not be obtained 
from the heat invariants ak{p,r) by radial expansions and linear combina- 
tions. That is, these curvature identities are not individually encoded into 
SIS{p). The explanation for this phenomenon is that they can be derived 
not from the original pre-functions defining the identities in SIS{p) but from 
functions defined by powering or multiplying the original heat invariant pre- 
functions. They usually define curvature identities lying in the complement 
of SIS{p). In the following discussions such equations are called powered 
curvature identities. 
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Although identities associated with f C^{u)dun or J C{u)H{u)dun do 
not show up in SIS{p), yet, on harmonic manifolds, they are determined by 
the spectra of local geodesic spheres. This statement is true for any iden- 
tity defined by J Ti{u)T2{u) . . . Tk{u)dun, where functions Ti{u) are density- 
Taylor coefficients obtained by expanding 9 regarding directions u. In- 
deed, on harmonic manifolds, the Taylor coefficients of f 6^{u,r)dun = 
(J 9p{u, r)dun)^ define spectrally determined identities, for any power /c G N, 
which are associated with integrals / Ti{u)T2{u) . . . Tk{u)dun- In order to set 
up a spectral data which contains also these powered density curvature iden- 
tities, system {afc(p, r)} should be extended by the functions f 6^{u,r)dun 
and extended spectral identity space, SIS^^{p), should be defined by this 
complete extended function-system. This space incorporates all identities 
which are determined by the spectra of local geodesic spheres of harmonic 
manifolds. 

It is interesting to see that how do these powered density functions define 
the powered curvature terms. We demonstrate this by considering the Taylor 
coefficient containing L in the expansion of 9. It can be proved that, by 
increasing power /c by 1, the contribution to the previous function is a linear 
combination of terms C^, CH, and L. Thus the powered identities associated 
with these pre-functions can be individually defined by powered density 
functions 9, 0^, and 9^. The AM&S paper uses normalized density for the 
expansions. In this paper we consider also the natural expansion ak{p,r) 
without normalizing the density 9 in the integral formula. By the above 
discussion, this exchange of normalized density to the non-normalized one 
defines new coefficients for the terms L,C^, and CH in the investigated 
Taylor coefficient associated with L. 

After these definitions we are able to explain why does the AM&S-proof 
breaks down after eliminating with the Lichnerowicz identity. The ultimate 
problem is that the curvature expression standing on the left side of this 
equation does not show up among the spectrally determined curvature ex- 
pressions. That is, this identity is not encoded neither in SIS{p) nor in 
SIS^^{p). This statement can be seen by proving that it is independent 
from both identity spaces. 

This observation is demonstrated, first, for the case satisfying 1 1 V-R| | (p) ^ 
0. Since it is independent from the higher order equations expressed in terms 
of higher order curvature terms, it is enough to see that it is linearly inde- 
pendent from the lower-order spectrally determined identities introduced 
so- far. In case of SIS{p), the problem reduces to the independence of (fT2]) 
from the system consisting (jlOp and the identity resulted by computations 
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([20 |) - (j23|) . The desired independence can immediately be decided by consid- 
ering only the main curvature expressions determined by terms R, R, and 
||Vi?|p. In case of SIS'^^{p), the data contains also other identities, but 
all new identities have main curvature expressions proportional to the main 
term of (fTUj) . Thus the independence is established also in this case. 

This independence from the whole SIS^^{p) means that the curvature 
expression on the left side of the Lichnerowicz identity can not be pro- 
duced by linear combinations of curvature expressions obtained from the 
expansions of ak{p,r) and 0^{p^r). It is obvious, that, the addition of a 
non-trivial constant-times of the Lichnerowicz identity to an identity be- 
longing to SIS'^^{p) changes it to one which does not belong SIS^^{p) any 
more. But exactly this happens during the considered elimination which 
can be described such that a certain linear combination of the Lichnerowicz 
identity and (fTOj) is added to the spectrally determined identity resulted in 
(I20])-(IM1)- Since ^ is in SIS''''{p), it stih keeps the identity subjected to 
elimination in SIS'^^{p). But the Lichnerowicz identity results an identity 
which is certainly not in SIS^^{p) any more. Therefore, this curvature iden- 
tity (expressed in terms of / L{u)dun, J C^{u)dun, f CH{u)dun, and ||Vi?||^ 
but not containing neither R nor R) can not be obtained by the radial ex- 
pansions of functions ak{p,r) and 6^{p,r). Since all terms but ||V-R|p can 
be obtained by such expansions, thus ||Vi2||^ is the only term which can 
not be computed by such expansions. This proves that it is not a spectrally 
determined quantity. 

In case of ||V-R|p = 0, the above investigated identities have only two 
main terms, R and R, thus the independence of the Lichnerowicz iden- 
tity from the other two identity is not insured by the above arguments. 
If it is independent, then one can imply in the same way that the elim- 
ination produces a spectrally undetermined identity. But such an iden- 
tity can not say anything about the spectral determinacy of local sym- 
metry. If it is in SIS^^{p), then terms R and R can be eliminated from 
the investigated identity. After this operation only terms associated with 
j C^{u)dunJ C{u)H{u)du n remain there. The only information provided 
by this step is that they are spectrally determined quantities. But this is a 
trivial information which does not imply the spectral determinacy of local 
symmetry, either. 

4.2 Rudimentary vs. refined SIS{p); Proper eliminations. 

Before examining the proof of the AM&S-statement in case of geodesic balls, 
we point out some of the rudimentary features of SIS{p). As it turnes out. 
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they can easily lead to misinterpretation of the meaning of the spectral 
detcrminacy of the curvature expressions defined by the given spectral data. 
Although the problems caused by these features can be fixed in case of 
geodesic spheres, but they are not so in case of geodesic balls. 

The primary problem about SIS{p) is that the curvature expressions 
obtained by evaluating integrals J Fc{u)dun, where Fc{u) is an c^'^-ordcr ex- 
pansion coefficient defined by radial expansion of a heat invariant ak{p,r), 
are considered without any reference to the degree c. It is obvious that this 
data can not be used for constructing functions depending on r. Such func- 
tions can be defined just by the complete expressions J Fc{u)dunr'^ involving 
r'^. This missing factor turnes the above elimination process to a very am- 
biguous operation. Indeed, if a term from J F(.{u)dun can be eliminated by 
a linear combination ^ • Qi f F^- {u)dun where none of the degrees Cj is equal 
to c, then that term is still present in J Fc{u)dunr'^ — Qi J Fc-{u)dunr'^'' ■ 
Furthermore, no linear combination A{p,r) = ^iPiakiiPif) exists whose 
coefficient of degree c is equal to the curvature expression resulting by the 
considered elimination. 

In order to explain this situation by a concrete example, suppose that 
density coefficient L appears in the coefficient (. . . , QL, . . . )n+i^""^^ of (n + 
l)-degree in the natural expansion of a heat invariant ak{p,r) what we 
wish to eliminate by the L appearing in the coefficient of (n -|- 5)-degree 
in the natural expansion of the density heat invariant J 9{p,r)dun- Then 
the coefficient (. . . , 0, . . . )„+i resulting by elimination can be decomposed 
as (. . . , 0, . . . = (. . . , QL, ■ ■ ■ — QL)n+i, meaning that it has ak{p, r)- 
and / r~^6{p, r)dM„-components. Since the latter one is not a heat invari- 
ant anymore, this elimination process involves new functions which are not 
listed among the elements of the spectral data and produces functions like 
L(u)(r"+^— r""*"^) which still contain L and which are not consistent with the 
requirement of radial constructibility of these functions. It is clear that all 
these problems are caused by eliminating terms appearing in coefficients of 
distinct radial degrees with each other. These arguments can not be brought 
up for eliminations performed by coefficients having equal degrees. Thus the 
only proper elimination in this situation is the elimination of QLn+i appear- 
ing in (. . . , QL, . . . )n+i by the Ln+i appearing in the corresponding coef- 
ficient defined by the natural expansion of f r^^9{p,r)dun = r~'^ J 9{p,r). 
Note that these coefficients are still spectrally determined by the spectra 
of local geodesic spheres, because, modulo certain multiplicative constants, 
just the degree defined by / 9{p,r)dun must be lowered by 4 in order to get 
the corresponding degree regarding J r~'^9{p,r)dun. Thus the proper treat- 
ment of these problems requires just the extension of the original spectral 
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data {ak{p,r)} to the data {ak{p,r), J{r''^0''{p,r)dun}- 

A precise establishment of this refined data, FineSIS{p), is as follows. 
In the first step the same type of expansion should be chosen regarding each 
heat invariant a^^p, r), which can be either the natural expansion or the nor- 
malized density expansion. The degree, c, associated with a spectral identity 
/ Fc{u)dun = {QR + . . . )c is indicated in the lower right corners. They are 
the coefficients of r'^'s appearing in the complete radial expansion formulas. 
For a fixed c, the SISc{p) is generated by all heat invariant coefficients of 
degree c. When this generator system is extended by coefficients, of degree c, 
defined for all functions J"(r^'^^*^(j9, r))du„}, then the corresponding extended 
coefficient space of degree c is denoted by SIS^-^"'- [p). In order to avoid 
operations among coefficients of different degrees, the total refined spectral 
identity spaces must be defined by the direct sums FineSIS{p) = (BcSISdp) 
resp. FineSIS^^'^{p) = ®cSIS^^'^{p). It is obvious that these spaces can 
be defined just for coefficients obtained from the expansions of ak{p,r) resp. 
J{r^'^e''{p, r))dun}. 

These refined spectral identity spaces put an end to the greatest imper- 
fection of the rudimentary S'/S(p)'s which do not give any indication about 
the radial functions from which the elements of the rudimentary SIS{p) are 
derived. In case of FineSIS{p) and FineSIS^^^ (p), however, any element 
J2i Qi f {u)d,Un of SISc{p) or SIS^^^{p) determines the expansion term 
Pc{r) = J2i Qi I Fc'\u)dunr'^ and there exist a function F{r) = '^^Fc{r) in 
the space spanned by functions {ak{p,r)} resp. {ak{p,r),J r^'^9^{p,r))dun} 
whose c*^ expansion term is Fc{r). Also the elimination process appears in a 
well defined clear form. If the above linear combination ^^Qi J Fc\u)dun 
eliminates a term then F{r) = '^Zc-^cir) ^^'^^ * globally defined function 
whose c*'* expansion coefficient is the curvature expression resulting by the 
elimination. Eliminations which can be described by such globally defined 
functions are called proper eliminations. Only such eliminations are accept- 
able to decide whether a term is spectrally determined by the given data. (It 
is obvious that all those curvature expressions are spectrally undetermined 
which can not be traced back to the functions {ak{p, r), f r'^'^9^{p, r)dun\ in 
this way.) Eliminations performed by the elements of SIS^^{p) are called 
rudimentary eliminations. Let it be pointed out again that, by the above 
arguments, these latter operations do not carry out the desired elimination, 
and, to the curvature expression, which are the results of formal rudimen- 
tary elimination, no globally defined function exists which satisfies the above 
properties. 
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4.3 Case of geodesic balls. 

A review of AM&S-computations. In |AM-S ] . the geodesic balls are 
considered regarding both the Dirichlet and Neumann conditions. These 
computations start out from the Branson-Gilkey formulas establishing the 
asymptotics of the Laplacian on a manifold with boundary. After adopting 
these generic formulas to harmonic manifolds, the Dirichlet coefficient a2{B) 
appears in terms of the shape operator au{r) = (S^)\T^xp(^ru)M (where u is 
the normal vector field on dB) in the following form: 

a^{B) = ^{ l^pf^(C\H)dvolB + J^^iPi'i-) + P3i.-))dvol9B), (28) 

where (C^ H) = 5(nC)^ - 2nC^ + ^n((n + 2)H - C^), (29) 
Pf (..) = 20nCTr{a) - 8CTr{a) + 16Tr{R^ o a), (30) 
Pf{..) = |(Tr(.))3 - ^lTr{a)Tr{a') + ^-^Tr{a% (31) 



The radial expansions of these functions are established by means of the 
following formulas: 

Tr(a) = (n- 1)- - -Cr - —Hr^ ^—Lr^ + O(r^), (32) 

^ > ^ V 3 45 15120 ^ ^' ^ ^ 

Tr{cj^) = (n - 1)1 - - ^^Hr^ + ^Lr^ + 0{r% (33) 

Tr{R, oa) = --- + {-^L + ^^TriUMy + . . . , (34) 

Tria^) = ■■■ + i^—L - —Tr(R'R')y + .... (35) 
^ ^ ^30240 96 ^ " ' 

In the last two formulas only the r^-coefhcients of the corresponding se- 
ries are indicated. Regarding the complete functions and P^, these 
coefficients can be written in the form 

[P,%{CH,L,Tr{R'M), [Ps%{C\CH, L,Tr{R'M), (36) 

where also the arguments (on which these coefficients depend) are indicated. 

Similar formulas are established also for a^{B). Although PfiC^, H) = 
P^{C^,H) and P^ = P^ , the third one: 

40 39 
P3^(..) = -^{Tr{a)f + 8Tr{a)Tr{a') + -Tr{a'). (37) 
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is linearly independent from Pf. Thus, also [P^]3{C^,CH,L,Tr{R'^R'J) 
is a new linear combination of the arguments. 

Normalized-density expansions and powered terms on balls. For 

the radial expansion of this local heat invariant, normalized density is used 
also in this part of the paper. It is computed just regarding the term defined 
on the boundary dB, in which case it means the expansion of 

•■ ^ ^^mW) I (-f"*-) +^3(..))<i-i*w). (38) 

In these discussions we need also the corresponding expansion of the volume 
function. By the previous formula, this means the expansion of 

^ ~^ — uJf \\ = ,/g / ^^ / / ^{Up, p)dupdp. (39) 
V0l{Sr{p)) V0l{Sr{p)) Jq Jsi(p) 

Note that this normalization concerns just the density defined at the end- 
point p = r of radial integration and it does not mean normalization regard- 
ing any p between and r. For explicit computations of the higher order 
powered terms appearing in the formulas one must consider also the so called 
higher order k-volumes defined by Volk{Bp{r)) = ^^^@'^{p, p)dvolBj,(r)i 
where G(r) = det A(cp(r)) = Aq + ^2?"^ + ^4?"^ + A^r^ + .... The exact 
formulas, regarding both the natural and normalized-density expansions of 
these functions, are described later in this section. 

Before going into the details let it be mentioned in advance that in case 
of geodesic spheres the spectrally determined functions Volk{Sp{r)) had to 
be extra added to the set {afc(p, r)} of local heat invariants, because, they 
do not belong to them originally, but, on harmonic manifolds, their spec- 
tral determinacy can be established by formulas such as / 6''{up,r)dun = 
{J 9{up,r)dUn)^ and / @''{up,r)dun = (/ @{up,r)dun)^. This spectral de- 
terminacy can not be established for functions Volk{Bp{r)), however. In- 
deed, for > 0, these functions are nothing to do with the heat invariants 
of local geodesic balls, furthermore, there does not exist a function-relation 
Volii.{Bp{r)) = Fi^{vol{Br{p)) expressing the k- volume Volk{Bp{r)) in terms 
of the standard volume vol[Bp{r)) = VolQ{Bp{r)), for all r. Thus neither 
the k-volumes nor the proper powered terms (obtained by using higher order 
k-volume-functions for eliminations) are determined by the spectra of local 
geodesic balls. Therefore this spectral data can not be extended by these 
spectrally undetermined functions. 

Rudimentary vs. proper eliminations on balls. After this detour we 
turn back to the AM&S-proof, where the desired spectral determinacy of 
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||Vi?|p is concluded just by a briefly explained argument. According to this 
explanation, just term L appearing in the above r^-coefiicients has to be 
eliminated by the L appearing in the corresponding coefficient of the radial 
expansion of the first heat invariant aQ{p,r) = vol{Br{p)) = {p,r). (The 
elimination of this L is really necessary because it involves rr(V„-R„V„i?„).) 
However, a deeper insight into the details shows that this briefiy explained 
step covers up numerous problems which discredit the AM&S-proof also in 
this case. One of them is indicated above. Even if this elimination were 
appropriate, there are still higher order spectrally undetermined powered 
terms left behind because of which term ||V-R|p can not be spectrally deter- 
mined. But in reality there are much more deeper reasons why this quantity 
is spectrally undetermined. 

First of all, the computations below show that the L appears with differ- 
ent degrees regarding the two functions ao{p,r) and a2{p, r), with respect to 
both type of expansions. The degree computed by a2ip,r) is always lower 
by 2 than the degree computed by ao (p, r) . This means that the proposed 
elimination is a typical rudimentary operation which leads to several con- 
tradictory situations. Remember that, in order to keep connection with 
the radial functions from which these coefficients are defined, this opera- 
tion must be well defined in terms of these functions. Since radial functions 
QiLr^ and Q2Lr'^^^ can not cancel out each other, the desired elimination 
of L docs not take place, furthermore, the coefficient defined by the for- 
mal rudimentary elimination docs not appear as a coefficient of a spectrally 
determined function, which, by definition, must be a linear combination of 
functions {ajt(p, r)}. In the top of these problems, also spectrally undeter- 
mined higher order powered terms are present in both functions. Since they 
can not cancel out each other, the proposed rudimentary operation is totally 
inappropriate for the desired elimination. 

Proper eliminations can be implemented by functions 



but this adjustment changes even vol{Bp{r)) into a spectrally undetermined 
function and it further worsens the spectral indeterminacy of the k- volumes 
of balls. Remember that, in case of geodesic spheres, the spectral determi- 
nacy of J r^6^{up,r)dun = J{9{up,r))^dun is due to the indicated equa- 
tion. Since is inside of an integration, such reconnection to the volume 
does not exist in case of geodesic balls. (By using multiple partial inte- 
grations, term p""^ can be retraced to p^~^. The complicated other terms 
produced by this process show then the impossibility of the reconnection of 




(40) 
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these adjusted functions to the volume.) In other words, all these functions 
are spectrally undetermined, together with all those terms which can be 
eliminated by these functions in a proper elimination process. This means 
that the L, which can not properly be eliminated from a2{p,r) by a spec- 
trally determined function, but just by a spectrally undetermined one, is 
spectrally determined in ao(p, r) whereas it is spectrally undetermined in 
'3^2 (Pi?^)- That is, the spectral determinacy of a term such as L is not an 
absolute but a relative concept. This property depends on the degree of the 
expansion term in which it shows up. This paradoxical statement points 
to the deep sitting ambiguity present in rudimentary spectral data defined 
for geodesic balls. It clearly explains also the error made in this part of 
|AM-S| . By the proposed elimination they tried to cancel out a spectrally 
undetermined quantity by a spectrally determined one. As it is explained 
above, this is a totally inappropriate operation from several point of view. 

Explicit expansion formulas. In this section those expansions are ex- 
plicitly computed which are mentioned in the above discussions. Expansion 
of vol{Br{p))/vol{Sr{p)) appears in the form: 

The details below show that this expansion does not provide L in a pure 
form. Instead, it is accompanied with higher order powered terms. In fact, 
a simple calculation shows that coefficients can be determined by the 
coefficients of expansions 

e(r) = det A(cp(r)) = Ao + ^2^^ + A^r^ + A^r^ + ..., (42) 

voliBrip)) = r"{-Ao + -^A2r^ + -^A^r' + -^A^r' + ...) (43) 
n n + 2 n + 4: n + b 

by the recursive formula 

= -^^^T^k - A2Vk-2 - A^Vk-A AuVq), (44) 



where k denotes even numbers. Thus (llOp really does not arise in a pure form 
in this expansion. Notice that, besides L, powered terms like and CH 
are also present in these local heat invariants. Thus the proposed rudimen- 
tary elimination contributes new proper powered terms to the ones already 
present there. Since the degrees defined by the considered two functions are 
different, these complications create new problems for the proposed rudi- 
mentary elimination. This phenomenon more emphatically underlines the 
inappropriateness of the proposed rudimentary elimination. 
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Also note that the expansion by normaUzed density represents these 
volumes Volk{Br{p)), where > 0, by series of the form r{Vko + Vk2r'^ + 
Vk^r^ + VkQr^ + . . . ) and such 0(r)-series appears just regarding the first 
function Pf {C^ , H)dvolB of a2{p,r). The other two terms define 0{^)- 
and 0(^)-series. The r^-component regarding the first function with re- 
spect to normalized-density expansion is a constant time of CPi{C'^,H), 
while this function is LPi{C'^ , H) regarding the r ''-component. 

The above arguments can easily be established also for natural expan- 
sions. Then function vol{Bp(r)) appears in the form described in (j43p . 
That is, no proper powered terms appear in this case which makes these 
computations more transparent than they were in case of normalized den- 
sity expansions. The natural expansion of the functions defined on the 
boundary dB can be established by expanding also 9 and multiplying this 
series with the AM&S-series. Since functions P^ and P^ have singular- 
ities of order O(^) and O(^) respectively, furthermore, f^P2 dvols in- 
volves integration regarding dr, this Taylor series can be written in the form 
Sfc(Qf(fe-4) + Q2(fc-2) + Q^k)^^'^^^ ■• where Qf^ is associated with the prod- 
uct of PP and e. That is: Qf = r"(Qfo + Qi2r^ + Q2 = r'^'^iQw + 
Qgr2 + .., and = r"-4(Qg) + Qg^^ _^ _ _ ^he above described r^- 
coefficients appear then in terms of Qi2{C^ ,CH), Q24{CH, L,Tr{R'^R'^)), 
and Qg{C^,CH,L,Tr{R'^R'J) where, comparing with those in [Pf]^, the 
arguments show up in completely new linear combinations. Indeed, the 
product of functions Pf^ with provides different contributions to the ar- 
guments of the individual terms labeled by i = 1,2,3. Unlike in case of 
normalized density expansions, no powered terms show up in the natural 
expansion formulas. But, because of the different degrees, the proposed 
rudimentary elimination still remains a completely inappropriate operation 
also regarding this expansion. 

Other explanations for the spectral indeterminacy of ||Vi?|p. In 
this section we describe two new interrelated constructions which also prove 
the investigated spectral indeterminacy. 

(1) In the previous constructions the required elimination of L and the 
other higher order density terms can properly be performed just by ad- 
justed volume functions. More precisely, this operation can be completed by 
the adjusted standard volume, Vol^'^^'^ ^ and two higher order k-volumes, 

Vol^'^^''' "* and Vol^'^'^'''' \ One can go in the other way around, by ad- 
justing a2(p, r) to the volume functions VoIq, Voli, and V0I2. This adjusted 
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function is then 

4'^*^'^''^')(5) = ^{ j^P?{C\H)dvolB + (45) 

+ / {T-'P^{..) + r^Pr,{..))dvol9B), (46) 

JdB 

where functions Pi, P2, and P3 are muhiphed with 1, r^, and respectively. 
This adjustment transforms the local heat invariant function to a spectrally 
undetermined function, which, like the volume function vol(B(r)), is of order 
0(r) in case of normalized density expansions and of order 0(r") in case of 
natural expansions. In each cases the L and the third order power terms 
defined on the boundary dB by the last two integrals involving P2 and P3 
appear with the same degree as L does regarding vol{B{r)). Also notice that 
the first integral P['{C^, H)dvolB defines term LPi{C'^ , H) regarding this 

degree. Term rr(V„i?„Vwi?u) present in this coefficient of 
can be obtained then by eliminating all these density terms by means of k- 
volume functions VoIq, Voli, V0I2, Vols, .... Because of term LPi(C^, H), 
this process can be completed by using more higher order volume functions 
than those used in the previous cases. Among these functions, which set 
includes also a^'^^'^'^ only VoIq is spectrally determined. Thus also 

this representation defines TriV uRu^ uRu) as a spectrally undetermined 
function. The same statement can be established by the adjusted functions 
^{adj,qx,q2r ,q3r where qi,q2,Q3 are arbitrary constants. 

(2) The above construction strongly suggests that the spectral deter- 
minacy of ||Vi?|p can be established by splitting a2{Bp{r)) into 3 diff'erent 
functions. The decomposition of a2{p,r) into three diff'erent parts is one of 
the most important characteristic features of these functions. The compo- 
nents in the coefficient {Q12 + Q24 + Q36) form a hnearly independent system 
of linear functions depending on the variables C'^,CH, L, and Tr{R'^R'^). 
System {L,Qi2,Q24,Q36} consists of for independent linear functions of 
these four variables, which appear as the corresponding Taylor coefficients 
of functions {ao{p, r), Qi{p, r), Q2{p, r), Qsip, r)}, respectively. By the given 
spectral data only 2 functions: 

{ao(p,r), ^(Qi(fc_4)(p) + Q2(k-2)(p) + Q3ik(p))r-"-^+'=} (47) 
k 

are determined by which no splitting of the second function can be estab- 
lished. It is the minimal data which provides the desired spectral determi- 
nacy of ||Vi?||2. 
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The same statements are true also for the spectral identities obtained by 
evaluating the integrals. That is, identity 

J {Qi2 + Q24 + Q36)du = iQi2iScal\Scal\\R\f)+ (48) 
+Q24iScal^,Scal\\Rf, - R,\\VR\\'^) + 

+Q36{Scal\Scal\\R\\'^, '^R - R,\\VR\f)) 

is the sum of 3 linearly independent identities which define a spectrally 
determined identity just in this linear combination. These component iden- 
tities together with spectral identity L{Scal^,Scal\\R\\'^,^R — R,\\VR\\'^) 
also determine a linearly independent system. The latter one is available 
now by the expansion of a^. The desired spectral determinacy of can 
be pointed out just by a drastically changed SIS{p) which, along with 
L{Scal^, ScalWRW^, ^R—R, ||Vi?|p), separately contains all above 3 spectral 
identities. That is, in case of geodesic balls, the split of the heat invariants 
into three parts provides a new spectral data by which the desired spectral 
determinacy can be established. Note that this is the minimal data by which 
the investigated problem can be positively answered. All attempts must fail 
unless they are able to establish this division of the heat invariant into 3 
separate invariants. But this split is certainly not encoded into the spectra 
of local geodesic balls, thus, all these attempts must fail altogether. 

Unified spectral data. Note that, regarding the two boundary condi- 
tions, the first two components in the above decompositions are the same, 
tliTis the difference eliminates both ones, leaving behind just — . 
Therefore, it seems to be possible to squeeze out some kind of positive an- 
swer regarding the spectral determinacy of ||V-R|p by eliminating L, Scal^, 
and 5*00/ 1 1 I p by using the two conditions together such that the unified 
data {a^{p,r),a^{p,r)} is extended by functions Volk{p,r). But also this 
elimination is just rudimentary which docs not work out by the very same 
reason explained earlier. In all of these cases proper elimination can be de- 
fined just by the refined spectral data FineSISB{p) = ®kSISB,k{p) which 
is the direct sum of SISB,k{py^ defined for radial degrees A; = 0, 1, 2, 3, ... . 
This definition excludes eliminations of terms with other ones which have 
different radial degrees. The positive conclusion can not be squeezed out 
even if the data is enlarged by the spectra of local geodesic spheres. In this 
case functions / 9^{p,r)dup exhibit contradictory properties regarding these 
two different type of spectra. Actually they are incompatible which is indi- 
cated by the distinct factors (47ri)~"'/^ and (47ri)~'^*'"*^'^)/^ appearing in the 
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asymptotic expansions of the partition functions, furthermore, the radial de- 
grees associated with the curvature terms are always even in case of geodesic 
spheres and odd in case of geodesic balls. Because of these differences, no 
proper eliminations of terms appearing in FineSISsip) can be defined by 
the terms appearing in FineSISg^{p). Thus the proper spectral data con- 
sidered for the unification of the ball-spectra with the sphere-spectra must 
be defined by the direct sum FineSISs{p) © FineSISB{p)- 

Summary 4.1. (A) Case of spheres. On harmonic manifolds the Lich- 
nerowicz identity Ci{p) is not encoded into the space, SIS^^{p), of spectrally 
determined identities defined for the spectra of local geodesic spheres. Since 
this space, defined by the radial expansion coefficients of functions ak{p,r) 
and j 0^{up,r)dup, is identical to the complete space of spectrally determined 
identities, the Lichnerowicz identity is not determined by the given spectral 
data, that is, it can not be computed by the spectra of local geodesic spheres. 
So is any identity T = S + QCi, where S S SIS^^{p) and Q ^ 0. In lAM-S^ . 
such spectrally undetermined identity is used to compute ||Vi?|p(p)|p. 

The greatest imperfection of SIS{p) is that the radial expansion coeffi- 
cients j Fc{u)dun, whose integral- evaluation provides the curvature expres- 
sions, are considered without any reference to the radial degree c, thus this 
data can not be used for constructing radial functions. (Such radial con- 
structions can be implemented just by functions J Fc{u)dunr'^ . Also note 
that, in case of geodesic spheres, this degree is even.) As it is pointed out 
above, eliminations using coefficients having different radial degrees result, 
on one hand radial functions from which the considered term is not elimi- 
nated, and, on the other hand, the curvature expression resulting by the for- 
mal elimination does not appear as an expansion coefficient of a function be- 
ing in Span{ak{p, r)} or Span{ak{p, r), J 9^{up, r)dunp}. In order to control 
this contradictory situation, refined spectral data FineSIS{p) = ®cSISc{p) 
resp. FineSIS^^{p) = ®cSIS^^{p), where SIS^p) resp. SIS^^{p) con- 
sists of coefficients having the same radial degree c should be introduced to 
the computations. This data is defined after choosing the same type of ex- 
pansion for each function from the set {ak{p,r)} resp. {f 6'^{up,r)dunp} . 
The direct sum applied in these definitions excludes the usage of coeffi- 
cients of different radial degrees during eliminations. Due to the identity 
f r'^'^9^{up,r)dup = r'^'^{f 9{up,r)dunp)^ valid on harmonic manifolds, a 
larger, spectrally determined data FineSIS^^{p) C FineSIS^^{p) can be 
defined by the elements of the larger function space 
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Although I110\) and the investigated r'^ -coefficients ofa2{p,r) appear with dif- 
ferent degrees, they can be brought together to have the same radial degree 
by the above degree- adjusting operations. Also note that the Lichnerowicz 
identity (which is not derived from a radial function l(p,r), thus no radial 
degree can be defined for it) can not be built into the space FineSIS{p). By 
the above arguments, it does not show up even in the rudimentary SIS{p). 
Therefore, the ultimate reason for ||Vi?|p(p)|p is not determined by the 
spectra of local geodesic spheres is that it can be computed just by the Lich- 
nerowicz identity. 

(B) Case of balls. In the ball case, the AM&S-paper proposes elimi- 
nation of L occurring in the r^ -coefficient of a2{Bp{r)) by the L occurring 
in the r"^ -coefficient of the spectrally determined function vol{Bp{r)) . Since 
these terms are defined for different radial degrees, by the above explana- 
tion, this rudimentary elimination can not be applied for these computa- 
tions. An other problem is that the considered -coefficient contains higher 
order powered terms which also must be eliminated. For this operation one 
must use also the proper k-volumes, Volk{p,r) = ^^-^ @pdBp{r) , defined 
for k > 1. Since these volumes are not determined by the standard volume 
vol{Bp{r)) = VolQ{p,r), they are not determined by the Dirichlet or Neu- 
mann spectra of balls either. Due to these facts, powered identities associated 
with J C'^{u)dun and f C{u)H{u)dun become spectrally undetermined. As a 
result, also ||Vi?||^ becomes spectrally undetermined. 

Proper elimination of these terms can be implemented just by using the 
adjusted volumes Vol'^-^''^ {p,f) = {r) P~'^^p^-^pi''') ' where k > 0. But 
these adjustments further worsen the spectral indeterminacy of the investi- 
gated terms. Even Vol'^'^'^'^ {p, r) becomes a spectrally undetermined quan- 
tity, meaning, that the spectral determinacy of L (regarding the spectra of lo- 
cal geodesic balls) is a relative concept. It depends on the degree of the radial 
expansion coefficient in which it shows up. In the considered r^ -coefficient 
it is a spectrally undetermined term and so are the other powered density- 
terms occurring there. This is an other proof of the spectral indeterminacy 
0/ ||Vi?|p. These arguments also show that, in case of geodesic balls, ex- 
tensions FineSIS^^{p) or FineSIS^^{p) provide spectrally undetermined 
quantities to the elements of FineSIS{p) . This is the most drastic differ- 
ence between the sphere- and ball-cases. 

Quantity ||Vi?[|^ can be determined just by appropriate extended data. 
Several of them are constructed above. All of them contain spectrally unde- 
termined identities which play the very same role than what the Lichnerowicz 
identity does when space SIS^^{p) of spectrally determined identities is ex- 
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tend by the Lichnerowicz identity in order to determine ||V-R|p. It is also 
pointed out above that this quantity is not determined even by any of the 
unified spectral data. 

5 Rudimentary uncertainty on Riemann manifolds. 

The above considerations exhibit strong physical contents. Namely, one can 
explain by the given spectral data that how docs the uncertainty princi- 
ple manifest itself on Riemannian manifolds. Here we consider only the 
rudimentary case. This idea works out on general Riemannian manifolds, 
however, by using the available formulas, we start to explain it on har- 
monic manifolds. Despite this analogy, it should be pointed out that this 
theory is pure geometric which can not be called physical, by any means. 
Although the constants and and relations associated below with spectrally 
undetermined identities (such as the Lichnerowicz identity) strongly remind 
the Planck constant and the corresponding Heisenbcrg relation, no physical 
interpretation for these identities exist at this point. This is why this theory 
is not associated with any kind of physics at this point. It has no apparent 
relation even to the exact physical models described in the last section. 

The constants and relations mentioned above appear on the scene when 
the distance of identities (such as the Lichnerowicz identity) from SIS{p) 
is computed. These computations are carried out in a Hilbert space in 
the following way. Consider all pre-functions F{u) defined in terms of 
arbitrary linear combinations of radial curvature terms appearing in the 
AM&S-expansions. In the expansions these terms appear just in particular 
combinations but in generalized situations all linear combinations are con- 
sidered. The integrals of functions F{u) define a linear space of equations 
endowed with the pre-Hilbert inner product defined by integrating products 
Fi[u)F2{u) by du. By topological closure, one defines an extended identity 
Hilbert-space EIS{p). Now let P{C) be the spectrally determined identity 
obtained by projecting the Lichnerowicz identity onto SIS{p). It can be 
obtained by integrating the projected function P{L{u)) by du. The identi- 
ties in SIS{p) are interpreted such that they can be clearly heard by the 
given spectral data while the other identities in EIS{p) which are in the 
complement of SIS{p) can be heard just together with noise. Such noisy 
identity is the Lichnerowicz identity, for which the noise wave is defined by 
NW{L{u)) = L{u) — P{L{u)). The magnitude of noise is defined by the L^- 
norm \\NW{L{u))\\^ = LiJ(||Vi2||'^, ...,RR, ...) which is equal to an expres- 
sion written in terms of the curvature terms appearing in the equations. Like 
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in the physical Heisenberg relations, these terms are typically products of 
those appearing in the primary equation MW[L). Note that these quantities 
are defined by the difference L{u) — P{L{u)), which makes the reminiscence 
with the Heisenberg relations even stronger. The Lichnerowicz-Planck con- 
stant, he, is defined by the minimum of values ||A?^VF(L('u))|p considered for 
all L{u) whose integration provides JC. Then the Lichnerowicz-Heisenberg 
relation is 

\\NW{Liu))f = LH{\\VR\\\\\VR\fR,\\VR\fR,RR, ...) > he, (49) 

which clearly describes the positive distance of C from the spectrally deter- 
mined identities. 

A more general version of these relations can be defined such that one 
replaces P{L{u)) by an arbitrary other function F{u) G Pre{SIS{p)), by 
which the noise wave, NWf{L{u)), is defined in the same way as by means of 
P{L{u)). Quantity ||A?^Wir(L(n))|p measures the magnitude of noise when 
the Lichnerowicz identity is listened by F{u). Formally this function looks 
the same as the previous one, but the noise is going to be bigger. That is, 
the left side of the general HL-relation ||A'"Wf(L(w))|P > he usually gets 
bigger in this situation. 

These objects can be defined on arbitrary Riemann manifolds for any 
identity. The spectral data can be defined by local geodesic spheres as well as 
balls. In order to avoid the above discussed problems, the AM&S-expansions 
must be established by the non-normalized density. Hilbert spaces SIS{p) 
and EIS{p) can be defined in the same way as they were on harmonic man- 
ifolds. It is also well known that Lichnerowicz established his identity for 
general Riemannian manifolds, where it appears in a much more complicated 
form, containing additional terms beyond those appearing on Einstein man- 
ifolds. Anyhow, one can prove also in these most general cases that this 
general identity is not encoded into the the spectra of local geodesic spheres 
or balls. More precisely we have: 

Summary 5.1. On a general Riemann manifold the Lichnerowicz identity, 
considered at a point p, is not determined by the spectra of local geodesic 
spheres or balls concentrated at p. That is, this identity is not encoded into 
the spectral identity space SIS{p) but it rather shows up as a spectrally unde- 
termined equation in EIS{p). This fact can be explained also such that pre- 
functions defining the Lichnerowicz identity do not show up in Pre{SIS{p)) 
hut in the larger pre- space ExPre{SI S{p)) = Pre{EIS{p)) where they have 
strictly positive distance from Pre{SIS{p)). The corresponding noise waves, 
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Plank constants, and Heisenberg relations can explicitly be computed also in 
these most general situations. 

This mathematical uncertainty theory matches the physical one, by all 
means. A curvature identity T defined by the integral J F{u)dun of a pre- 
function is always a true identity. This logical value is independent from its 
relation to the spectral identity space SIS{p). The above quantities measure 
that in what extend can this well determined identity be measured ( or, recov- 
ered) by the given spectral data. It is also clear that spectrally undetermined 
identities come from two sources. Either by powering spectrally determined 
identities, or by considering such identities like Lichnerowicz's which are 
nothing to do with pre-functions F{u) G Pre{EIS{p)). 

Several alternative versions of the above concepts can be introduced as 
follows. For establishing the first one, consider a geodesic sphere or ball 
of radius R{p) about the center p. In the following definition the data can 
either be the spectra of local geodesic spheres or balls concentrated at p. If 
the Taylor expansion a^ip, r) = Yl'^o h^ks{p)f"'^ of the local heat invariants 
is convergent for r = R{p), then, for any fixed k, heat invariant ak{p,R) = 
^,^ks{p)R^ (of the sphere or ball of radius R{p)) can be considered 
as a curvature identity, Ak{R{p)), defined as an infinite sum (series) of 
identities -fr^fcs € SIS{p) determined for the coefficients of this expansion. 
By this construction, one can associate the identity space SIS{R{p)) = 
Yl^=o -^kiRip)) to the spectrum of any fixed geodesic sphere or ball, which 
is just a very thin subspace of the total space SIS{p). For a field of geodesic 
spheres or balls, having radius R{p) at a point p, the SIS{R{p)) defines 
a Hilbert space bundle over the manifold, for which the Planck constants 
and Heisenberg relations can be defined in the same way as for SIS{p). If 
an identity is not determined by SIS{p), it is even so regarding SIS(R{p)). 
This construction shows that the information encoded into the spectral data 
is considerably weakened if, instead of all local geodesic spheres or balls 
concentrated at p, it is provided only by a single sphere or ball having 
center p. 

Other concepts analogous to the AM&S local heat invariants are the so 
called remote local heat invariants, pkip, t), defined at a point p of a compact 
Riemann manifold by means of the cut locus Cut{p). For a unit vector Up G 
Tp(M), let q{up) G Cut{p) be the point such that the geodesic starting out 
from p into the direction of Up intersects Cut{p) at q{up). Parameterization, 
r, on this geodesic, Cpqir), is chosen such that Cpg(O) = p, Cpg(l) = q{u), 
furthermore, the speed vector Cpq{T) has constant length equal to the arc- 
length of Cpq. This parameter is called cut-locus-radius function. For a fixed 
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< r < 1, cut-locus-sphere CuSp{T) of center p and radius r is defined by 
consisting of points which have parameter r on each of the above geodesies. 
Cut-locus-balls, CuBp{T) are defined by points having parameters < r. The 
remote local heat invariants are defined by the cut-locus-radial expansion 
(i. e., r-expansion) of the heat invariants afc(CwS'p(r)) resp. ak{CuBp{T)) 
defined on cut-locus-spheres resp. cut-locus-balls. Spectral identity spaces 
SIS{Cut{p)) and SIS{Cut{T{p))), where r denotes fixed cut-locus-radius, 
regarding cut-locus-spheres resp. balls can be introduced in the same way 
as for geodesic spheres resp. balls. 

The above constructions allow to introduce spectral identity spaces also 
for compact sub-domains, D, with boundaries, dD, of Riemannian manifolds 
M". For the sake of simplicity suppose that dD is diffeomorphic to the 
Euclidean unit sphere S^~^, furthermore, let p G D be a point such that 
for any q € dD there exists a unique geodesic Cpq C D joining the two 
points. On such a geodesic, one can define the same parameterization, Tp, 
satisfying the very same properties described above. By replacing cut locus 
Cut{p) with boundary dD in the above constructions, one can define spectral 
identity space SIS{p, dD) (where dD appears as level set Tp = 1) regarding 
the spectrum of Laplacian defined on dD or any of the spectra defined by 
particular boundary conditions on D. By continuous movement of p, one 
can define this identity space together with the associated Planck constants 
and Heisenberg relations at any point lying in the interior of D. Both the 
spectrally determined and undetermined identities can precisely be described 
in terms of the Hilbert space bundles established in this construction. 

6 Global vs. local spectral investigations. 

On compact Riemann manifolds, M, the integrals of the infinitesimal heat 
invariants define the so called averaged infinitesimal heat invariants. The 
exploration of relations between the global spectra and the averaged in- 
finitesimal heat invariants seems to be an interesting question. In this paper 
we establish just a single theorem concerning the relation between the global 
spectra and the volume of small geodesic balls and spheres. Next we prove 
that, contrary to the volume of the whole compact manifolds, these local 
averaged volumes are not determined by the spectrum of M. 

By ([2]) , the fourth coefficient appearing in the radial expansion of ^^^j,^ is 
a linear combination of terms TrR^^ and TrR'^^ = Ricc{up, Up)\up\up, where 

1 Up I Up denotes two covariant differentiations regarding Up. Integration of 
TrR'^ regarding dup at a fixed points provides linear combinations of terms 
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like A{Scal)(p) and {Ricc'^l^)^^, whose integrals vanish on the whole compact 
manifold, by the Stokes theorem. The same integral applied to the first term 
provides a constant time of ||i?zcc|p + (3/2)||i?|p. 

To establish the statement we consider also the one parametric families 
{G, gx[t) ) of 0-isospectral manifolds constructed by Schiith [Sj on compact 
Lie groups G such as: 

SO{m) xT"^ {m> 5), Spin{m) x {m > 5), SU{m) x T'^ {m > 3), (50) 
SO{m) {m > 9), Spin{m) (m > 9), SUim) (m > 6), SO{8), Spin{8). (51) 

She proved that distinct members in a one-parametric family have distinct 
spectra on 1-forms by showing that they have Ricci tensors of different norm. 
We use the very same computational technique (she applied to prove this 
latter statement) for establishing our statement. 

Consider a2{G, g^^t)) for ^ family, which, due to the isospectrality, is 
constant regarding t. Since Seal is constant on G, furthermore, vol{G) 
and f^Scal are spectrally determined, also the first term, J^Scaf', in 
a2{G, 5A(t)) is constant regarding t. It follows then that the averaged volumes 
of small geodesic balls and spheres can not be constant, otherwise the above 
fourth coefficient in the expansion of the density together with a2{G,g\(t)) 
define a non-degenerated system of linear equations which would determine 
constant values both for ||fficc|p(t) and ||i?|p(t). This idea works out for 
all invariants TrA^^^y By summing up we have 

Theorem 6.1. The averaged volumes of small geodesic balls and spheres 
are generically not determined by the spectra of compact Riemann mani- 
folds. This spectral indeterminacy is exhibited by the one parametric fami- 
lies {G,gx(t)) constructed by Schiith JSj^ on the compact Lie groups G listed 
in /i50\) and Ii51\) . The fourth coefficient in the power series expansion of the 
volumes of small geodesic balls resp. spheres depends on parameter t. This 
coefficient is expressed in terms of \\Ricc\^ and other lower order curva- 
ture invariants. Schiith originally proved that quantity ||i?zcc|p is different 
on the distinct members of a family, therefore, she concluded, they are not 
isospectral on 1-forms. 

This statement can be established, by the same proof, for TrA'^-volumes 
of geodesic balls and spheres where the integration can be defined either 
with 9du or just with du. This statement is highly expectable also on the 
Gordon- Wilson examples |GW] of isospectral manifolds having different lo- 
cal geometries as well as on those derived from these examples. Since these 
manifolds have non-constant scalar curvature in general, this version of the 
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theorem can not be established by using only the fourth coefficient of the 
above power series expansion. 

7 Isospectralities and physical symmetries. 

Both my isospectrality examples and the Gordon- Wilson examples |GW) 
(together with those derived from the GW-constructions, see a survey on 
them in [S]) arise from 2-step nilpotent Lie groups. Yet, in terms of spectral 
stability of small geodesic balls and spheres these two types of examples are 
the polar opposites of each other. One of the idiosyncrasies of the GW- 
examples and their relatives is the wide range of changing invariants during 
the continuous isospectral deformations introduced in their constructions. 
Whereas, these quantities are not changing during the discrete isospectral 
deformations introduced in my constructions. In my isospectrality families 
the members share even the same volume function 0(r), therefore, also the 
volumes of geodesies spheres resp. balls having the same radius must be 
identical. Thus the above theorem is a clear demonstration of these argu- 
ments. Later, after introducing the technical definitions, it is more clearly 
pointed out that what causes these significant differences between these two 
different type of constructions. 

The isospectrality exhibited by my examples has a deep physical mean- 
ing. Namely, it is equivalent to the C-symmetry known in quantum theory. 
This physical connection is completely established in |Sz8] . Actually, the 
bulk of this whole section is going to be an outline of some of the results 
developed in |Sz8j . Readers interested in details may consult with this pa- 
per, however, the following review tends to be self-contained as much as it 
is possible. 

Before going into technical details, we briefly describe this physical con- 
nection by non-technical words. As it is pointed out below, on the Rie- 
mann manifolds used for these constructions the Laplacian is nothing but 
the quantum Hamilton operator of particle systems described in elementary 
particle physics. The action of the intertwining operator constructed for es- 
tablishing the isospectrality can be interpreted such that it exchanges some 
of the particles for their anti particles. Thus the mathematically established 
isospectrality can be paraphrased such that the spectrum of an elementary 
particle system, that is, the possible energy levels on which the system can 
exist, is not changing if some of the particles are exchanged for their anti 
particles. 

This paraphrase is a clear manifestation of the C-symmetry principle 
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introduced in physics in the following general form: "The laws are the same 
for particles and antiparticles." It should be emphasized, however, that this 
general principle will be pointed out only regarding the spectrum of the 
Hamilton operator. In other words, the isospectralities exhibited by my 
examples are clear manifestation of the spectral C-symmetry principle. 

A preliminary description of the endomorphisms exchanging the parti- 
cles for their antiparticles is as follows. The two-step nilpotent Lie groups 
used for these constructions are defined on the (X,Z)-space, M!' xM}' = X x Z 
by a linear space, Eshewi''^)^ of skew endomorphisms acting on the X-space. 
The members of an isospectrality family are defined on the same (X,Z)- 
space by a corresponding family of endomorphism spaces where for any two 
members, Eg^^^^X) and E'gf^^^{X)., there exists an involutive orthogonal 
transformation a : X ^ X commuting with all these endomorphisms and 
the endomorphisms belonging to E'^^^^{X) can be obtained by composing 
the elements of Eskew{X) by cr, that is, E'^k^^{X) = a o Eskewi^)- These 
conditions imply the existence of an orthogonal direct sum 
where both subspaces are invariant under the actions of all endomorphisms 
and the first component is fixed under the action of cr, while it is —id^(b) on 
the orthogonal complement X^^\ It is explained later that these endomor- 
phisms serve as angular momenta for particles orbiting in X^"''^ resp. X^^\ 
These endomorphisms reverse the angular momentum for particles living on 
;f(^) and preserve it for those orbiting in X^"'\ The operators intertwining 
the Laplacians arc defined by means of these exchange-endomorphisms. Let 
it be mentioned yet that these constructions allow just discrete isospectral 
deformations. 

By using the same two-step nilpotent Lie group in two different ways, 
these examples are constructed on two different type of manifolds. The first 
ones are torus bundles defined over the X-space by factoring the center, Z, 
by a Z-lattice F^, and the others are ball resp. sphere bundles over the 
X-space defined by considering appropriate Z-balls resp. Z-sphcres in Z. 
For the sake of simplicity, next wc describe only the ball-bundle cases. The 
spectral investigations must be established on the function spaces defined 
on these manifolds. Actually both function spaces can be considered on 
the same non-factorized manifold. Namely, in the first case, it is the space 
of F^-periodic functions, while in the second case it is the function space 
satisfying an arbitrarily prescribed boundary condition on the boundary 
manifold which is a Z-sphere bundle. Both function space can be written 
up by appropriate Z-Fourier transforms. In the first case it is the discrete 
Z-Fourier transform defined by the Z-lattice. The big advantage of this 
transform is that it separates the X- and Z-variables from each other. The 
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other Z-Fourier transform introduced for investigating Z-ball bundles is the 
so called twisted Z-Fourier transform. This name indicates that this is a 
much more comphcated version of the Z-Fourier transforms where the X- and 
Z-variables can not be separated from each other. This transform involves 
mixed, so called twisting functions which establish the fundamental bonds 
between the X- and Z-spaces. 

One of the most important observation regarding these two exact mathe- 
matical models is that with the help of them the electromagnetic, the weak-, 
and the strong-nuclear forces can be established within a unified frame- 
work. The main unifying principle is that these forces can be described by 
the eigenfunctions of the very same Laplacian such that the distinct forces 
emerge on distinct invariant subspaces of this common quantum operator. 
This very same Laplacian is the Laplacian acting on the Lie group and the 
two distinct function spaces are the F^-periodic function space resp. those 
defined for Z-ball bundles. After its action on the Fourier integral formula, 
the very same Laplacian appears in drastically different ways behind the 
integral sign. In the first case, it turns into a Ginsburg-Landau-Zeeman 
operator of charged particles orbiting in complex planes in magnetic fields 
perpendicular to the planes, while, in the second case, it appears as a sum 
of a scalar operator which can be identified as the quantum Hamilton oper- 
ator of electroweak interactions and a new type of spin operator which can 
be identified as the quantum Hamilton operator of strong force interaction 
keeping the nucleus together. The particles attached to the two function 
spaces are distinguished by calling them particles having no interior resp. 
those having interior. 

The Hamilton operators are defined in this theory on the nilpotent 
groups. The corresponding wave and Schrodinger operators emerge in the 
Laplacians of the static resp. solvable extensions of these nilpotent groups. 
The latter manifolds are endowed with a natural invariant indefinite metric 
of Lorentz signature. Thus, these new exact mathematical models provide 
a relativistic theory for elementary particles. The above discussed functions 
defined by Z-Fourier transforms appear in the explicit solutions of these 
wave operators. These wave functions strongly remind those introduced by 
de Broglie in classical wave mechanics (cf. [P], volume 5). Actually, they 
are the only appropriate adaptations of the original de Broglie wave func- 
tions to the new mathematical models. They carry over also the original 
de Broglie theory to understand the much more complicated physical situ- 
ation inherent in the new models. This theory establishes infinitely many 
non-equivalent models for which even classification is possible. The parti- 
cle systems attached to them behave exactly like those introduced by the 
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familiar standard model of elementary particle physics. 

7.1 Technicalities on 2-step nilpotent Lie groups. 

A 2-step nilpotent metric Lie algebra, {M, {,)}, is defined on a real vector 
space endowed with a positive definite inner product. The name indicates 
that the center, Z, can be reached by a single application of the Lie bracket, 
thus its second application always results zero. The orthogonal comple- 
ment of the center is denoted by X. The Lie bracket operates among these 
subspaces according to the following formulas: 

[M,M] = Z , [M,Z]=0 , M = X®Z = R'' xRK (52) 

Spaces Z and X are called Z- and X-space, respectively. 

Up to isomorphisms, such a Lie algebra is uniquely determined by the 
linear space, Jz, of skew endomorphisms Jz '■ X ^ X defined for Z-vectors 
Z £ Z hy the formula 

{[X,Y],Z) = {Jz{X),Y),yZ GZ. (53) 

This statement means that for an orthogonal direct sum, Af = X (B Z = 
R'^ X M', a non-degenerated linear map, J : ^ — > Eskew{X) , Z ^ Jz, from 
the Z-space into the space of skew endomorphisms acting on the X-space, 
defines a 2-step nilpotent metric Lie algebra on TV by (I53p . Furthermore, 
an other non-degenerated linear map J having the same range Jz = Jz as 
J defines a Lie algebra which is isomorphic to the previous one. If isometric 
isomorphism is required, then o J must be an orthogonal transformation 
on the Z-space. 

With the help of these technical definitions one can more clearly explain 
why do the GW-type constructions produce completely different examples 
from those I constructed by the above exchange-endomorphism a. In the 
GW-case the isospectral deformation is defined such that, for any fixed Z, 
the spectrum of the bilinear map (Jz(Xi), Jz{X2)) (which defines a unique 
linear map on the Euclidean X-space) is not changing, but, the spectrum of 
(J^^ (X), (^)) (defined, now, on the Euclidean Z-space for an arbitrary 
fixed X) is changing. Other characteristic features are that these maps have 
different eigenvalues, furthermore, the Z-space has dimension 2. The chang- 
ing spectrum of {Jzi{X), Jz^iX)) gives rise to the wide range of changing 
invariants (including also spectral invariants) defined for Z-balls, Z-spheres, 
geodesic balls, and geodesic spheres. This is the ultimate reason why the 
GW isospectrality examples are established just on Z-torus bundles but not 
on Z-ball and Z-sphere bundles, where they are actually non-isospectral. 
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Whereas, the spectra both of {Jzi^i), Jz{X2)) and ( J^^ (X), (^)) 
are not changing during the discrete deformations I introduced for the con- 
structions. This is why my isospectrahty examples hve both on Z-torus 
bundles and Z-ball resp. Z-sphere bundles. Due to these differences, the 
following formulas established in my case resist any attempt to define them 
on the GW-examples. In this paper we describe only such cases where the 
spectra of the above maps consist only the same values. Thus difficulty 
arises already in the first step when one tries to extend them to the GW- 
type constructions, where these spectra consist of different eigenvalues. But 
the ultimate reason for this extension breaks down is that the spectrum of 
{J z-i{X) , J Z2{X)) is changing during the continuous deformations they in- 
troduced for their constructions. These difficulties strongly indicate that the 
distinguishing characteristics exhibited by these two different cases can be 
explored by mutually distinct methods. This field can not be explored by a 
single method with the help of which one would be able to control both the 
changing and the non-changing invariants exhibited by these two different 
type of constructions. 

Important particular 2-step nilpotent Lie groups are the Heisenberg-type 
groups [K] defined by endomorphism spaces Jz satisfying the Clifford con- 
dition = —z'^id, where z = \Z\ denotes the length of Z-vectors. These 
groups are attached to Clifford modules (representations of Clifford alge- 
bras). The well known classification of these modules provides classification 
also for the Heisenberg-type groups. According to this classification, the 
X-space of a H-type group is an (a + 6)-times Cartesian product of a smaller 
space y = M"', which is endowed with an ^-dimensional endomorphism 
space jz such that the endomorphisms acting on X = y"" + are defined 
by = jz X • • • X jz X —jz X • • • X —jz- Dimension ni depends only on /. 

By means of the exponential map, also the group can be considered to 
be defined on M. That is, a point is denoted by {X, Z) also on the group. 
Metric tensor, g, is defined by the left invariant extension of (, ) onto the 
group N. 

Although most of the results described below extend to general 2-step 
metric nilpotent Lie groups, next only H-type groups will be considered. On 
these groups, the Laplacian appears as follows: 

1 

A = Ax + {l + -\X\^)Az + ^daDa; (54) 

a=l 

where {bq} is an orthonormal basis on the center (Z-space) and Da* de- 
notes directional derivatives along the vector fields X Ja{X) = Je^{X), 
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furthermore, x = |X| denotes the length of X- vectors. 

In the next sections, it is described in technical terms that how this 
operator manifests itself as the Hamilton operator of different elementary 
particle systems. A non-technical preview of this physical interpretation is 
as follows. As it is indicated above, the main idea is that these systems are 
attached to various invariant function subspaces of A which can be divided 
into two major classes. The first class is defined by Z-torus bundles (alias Z- 
crystals), where the corresponding function space consists of functions which 
are periodic regarding the Z-lattice Tz defining the Z-torus bundle. The par- 
ticles attached to such a function spaces are considered to be point particles 
having no interior. They can show up at the lattice points of Tz, where 
the angular momentum is defined for them by means of Jz- As it will turn 
out, the A appears on this function space as the Ginsburg-Landau-Zeeman 
operator of charged particles. Thus the natural interpretation for these 
systems is that they consist of electrons, positrons, and electron-positron- 
neutrinos. Contrary to these cases, the other type of particles, emerging on 
function spaces defined on Z-ball bundles by Dirichlet or Neumann bound- 
ary conditions, do have interior which space is represented by the interior of 
the Z-balls. The X-space represents always the exterior word. The Lapla- 
cian appears on these Z-ball-bundles as the sum of Hamilton operators of 
electro- weak and strong-force interactions. A technical description of these 
enormous differences between Z-torus and Z-ball bundles is as follows. 

7.2 Z-crystals modelling Ginsburg-Landau-Zeeman operators. 

The Z-torus bundles are defined by factorizing, T\H, a two-step nilpotent 
group H hy a Z-lattice, F = {Z-y}. The name indicates that this lattice is 
defined only on Z and not on the whole {X, Z)-space. Such a factorization 
defines a Z-torus bundle over the X-space. The natural Z-Fourier decompo- 
sition, := J^-y of the function space is defined such that subspace 
is spanned by functions of the form 



Note that each Wj is invariant under the action of A, more precisely, we 
have: 



Z) = 'ilj{X)e 



,2Tri{Z^,Z) 



(55) 



A'^^{X,Z) = {<i^xp){X)e^'''^^-"^\ where 
<i-y = Ax + 2iTiD^ • -47r2|Z^p(l + 



(56) 
(57) 
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In terms of parameter /x = 'k\Z\^^ this operator can be written in the form 
<l^ = Ax + 2\D^ • — ;U^|Xp — Although it is defined in terms of the 
X-variable, this operator is not a sub-Laplacian which is resulted by a sub- 
mersion. It should be considered as restriction of the total Laplacian onto 
the invariant subspace W^. Actually, the Z-space is represented by the con- 
stant [1 and operator A characteristic feature of this restricted operator 
is that it involves only a single endomorphism, Jz^ ■ 

On a Z-crystal, BrxT\ for given boundary condition Af'{B?)+Bf{B?) = 
0, the eigenfunctions of <l^ can be represented in terms of n^'^-order com- 
plex valued spherical harmonics h("'''")(X) in the form /((X, X))h("'''") (X), 
where the radial function / is an eigenfunction of the ordinary differential 
operator 

(<>^,t/)W =4t7''W + (2fc + 4n)/'(tV(2m/x + + (58) 

Degree m is defined such that h'^"''™)(X) is simultaneously eigenfunction also 
of with eigenvalue m/i. 

Polynomials h('''™)(X) can be constructed as follows. Consider a com- 
plex basis B = {Qi, • • • , Qfc/2} regarding the complex structure J^^„ and an 
n*''-order polynomial \[zf{Z^u,X)'zl\Z^u,X), where Zi{Z^u,X) = {Qi + 
^JZf^uiQi):^) and YliPi + Qi) = n. Since 2;j(Z^„,X) is an eigenfunction of 
i-D^» with eigenvalue — ;U, the whole polynomial is also an eigenfunction with 
eigenvalue m/i, where Yli~Pi + = However, this homogeneous poly- 
nomials are not harmonic regarding Ax and their restrictions do not define 
spherical harmonics on the unit X-sphere. In order to have the spherical 
harmonics, these spherical but non-harmonic functions should be projected 
into the space of n*^-order spherical harmonics. These projections, II^\ 
are extensively described in my papers [Sz4j - [Sz8] (see for instance Section 
6.4 in |Sz8] ). They can be represented as uniquely determined polynomi- 
als of the spherical Laplacian A5. Since operator iD^» commutes with 
this Laplacian, projection n^^((]^ 2;f'z?'))(Z^„, X) really provides a desired 
polynomial h("''")(X). Since these projections are onto whose kernels are 
formed exactly by the lower order polynomials, all these polynomials can be 
obtained by this construction. 

In the 2D-case, operator <^ can be transformed to the the well known 
Ginsburg-Landau-Zeeman operator \Bo\ \LL\ IP] 

-i^\x,y)-i^D,*+—^{x^ + y^) (59) 
Am ^ ' Zmci omc^ 
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of a charged particle orbiting in the (x, y)-plane in constant magnetic field 
directed toward the z-axis by choosing /i = eB/2hc and multiplying the 
whole operator with —h"^ /2m. For k-dimensional X-spaces, number k = k/2 
means the number of particles, and, endomorphisms and —jz in the above 
formulas are attached to systems electrons resp. positrons. More precisely, 
by the classification of H-type groups, these endomorphisms are acting on 
the irreducible subspaces M"' and the system is interpreted such that there 
are n/ /2 particles of the same charge orbiting on complex planes determined 
by the complex structures jz^ resp. —jzu i^i constant magnetic fields. 

Note that this operator contains also an extra constant term 4^^, which, 
after establishing the corresponding wave operators, can clearly be explained 
as the total energy of neutrinos accompanying the electron-positron sys- 
tem |Sz8j . (This energy term is neglected in the original Ginsburg-Landau- 
Zeeman Hamiltonian.) Thus these mathematical models really represent 
systems of charged particles and antiparticles accompanied with electron- 
positron-neutrinos known in elementary particle physics jV l IWll IW2) . Let 
it be mentioned yet that operator is associated with the magnetic dipole 
resp. angular momentum operators of classical quantum theory. They are 
defined for the lattice points separately. Thus, while wandering on the lat- 
tice points, these point-like particles receive their angular momenta at the 
lattice points where they stay on. 

Finally, let the isospectrality question be clarified. An isospectrality 
family is defined by Heisenberg type groups Hi""'^^ having the same / and 
(a + 5), that is, they share the same X-space <V = y"" + and Z-space 
Z. They differ from each other just by the decomposition of the X-space 
and the action of Jz on these components. When }{^°'~^^'^^ ig compared 

with Hi "''^\ then the exchange endomorphism a described above is defined 
by id/ya resp. —id^yb on the components of the above decomposition. On 
^(b) = y'', this exchange endomorphism defines the angular momentum 
of H^""'^^ by the negative of the angular momentum defined for }{^"'~^^'^\ 
One can interpret this as switching the sign of charge of the particles on 
this component. Now consider an X-ball around the origin and restrict 
the torus bundle onto this ball. For both bundles, the normal vector at a 
boundary point is the the radial unit X-vector. For both metrics consider 
the same complex basis B and the same polynomials in terms of the complex 
structures defined for these two metrics. Since both define the same radial 
Laplacian with the same boundary conditions, these two metrics on the 
considered sub-bundles must be isospectral regarding any of the boundary 
conditions. This isospectrality can be seen also by observing that, for any 
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fixed Z^, there exist an orthogonal transformation on the X-space which 

conjugates J^""''^'''^ *° Jz^'^\ therefore, it intertwines the Laplacians of the 
two metrics term by term along with the boundary conditions. 

If basis is chosen by picking the vectors always from subspaces 
y = M"'^'^ this isospectrality is a clear manifestation of the spectral C- 
symmetry. For a general basis, however, it exhibits also internal symmetries. 
In elementary particle physics, this name was chosen to indicate that one 
think about internal symmetries as having to do with the intrinsic nature 
of the particles, rather than their position or motion. You can think of 
each particle as carrying a little dial, with pointer that points in directions 
marked "electron" or "neutrino" or "photon" or "W" or anywhere in be- 
tween. The internal symmetry says that the laws of nature take the same 
form if we rotate the markings on these dials in certain ways. If a basis 
vector Qi = Q^f^ + Qf^ is lying neither in 3^*-"^ nor in y^'^\ the exchange 

endomorphism effects only Qf \ that is, the particle antiparticle exchange 
is just partial and not complete. By this explanation, this isospectrality is 
a clear manifestation of the spectral C/I-symmetry. 

7.3 Extended particles occupying Z-ball bundles. 

Contrary to the above Z-crystal cases, the Laplacian on Z-ball and Z-sphere 
bundles can be identified with the Hamilton operators of particles to which 
interior can be attributed. In elementary particle physics such particles 
occur in the nucleus where the constituents are held together by the elec- 
troweak and strong forces. These forces are exhibited by the idiosyncratic 
appearance of the Laplacian on these bundles. Namely, it decomposes into 
a scalar operator and a non-standard spin operator (called also roulette op- 
erator) where the scalar operator represents the electroweak force while the 
roulette operator corresponds to the strong force. 

The ball X ball- and ball x sphere- type manifolds used to these investiga- 
tions emerged first in the spectral constructions performed in [Sz4] - [Sz6] . 
These manifolds are defined by appropriate smooth fields of Z-balls resp. Z- 
spheres of radius iJzd-'^l) over the points of a fixed X-ball Bx of radius Rx- 
Note that radius i?z(|X|) depends just on the length, \X\, of vectors X S Bx 
over which the Z-balls resp. Z-spheres are considered. The boundaries of 
these manifolds are the so called spherexball- resp. sphere x sphere- type 
manifolds. Comparing with the Z-crystals, the difference between the two 
type of bundles is that one considers Z-balls resp. Z-spheres instead of the 
Z-tori used in the previous construction. In the isospectrality investigations 
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the compact domains corresponding to i?^ < oo play the primary interest. 
In physics, however, the non-compact bundles corresponding to Rx = oo 
(that is, which are defined over the whole X-space) become the most impor- 
tant cases. In the following considerations both cases will be investigated. 
The details will be provided in this paper just for Z-ball bundles and not for 
Z-spheres bundles. 

The main difference between the Z-crystals and Z-ball bundles is that the 
computations in the latter case can not be reduced to a single endomorphism 
but they must be established for the complete operator M = ^ daDa*- This 
operator includes the angular momentum endomorphisms Jz with respect 
to any Z-direction. This complication gives rise to a much more complex 
mathematical and physical situation where both the exterior and the interior 
life of particle systems exhibit them-self on a full scale. 

The above argument implies that this mathematical situation can not 
be described by the discrete Z-Fourier transform applied on Z-crystals. In 
this case one considers a fixed complex basis B together with the complex 
coordinate systems {zi{Ku, X)} defined, on the X-space, by the complex 
structures Ji^„. Then, the Z-Fourier transform is defined on the whole cen- 
ter, Z = M', by 



It should be pointed out that a fixed B can serve as a complex basis only 
for almost every Ku, which form an everywhere dense open subset on the 
unit Z-sphere. However, the polynomials are well defined analytic functions 
even at those K^s with respect to which the B is not a complex basis. This 
formula is well defined if for any fixed \X\ function j4(|X|,if) is of class 
regarding the i^-variable. 

This so called twisted Z-Fourier transform does not separate but rather 
binds the Z-space and the X-space together. This strong bond is established 
by the polynomials Y\ -^f*^?* which depend both on the X- and K-variables. 
It is proved in Section 6.2 of [Sz8] that the complex valued functions consid- 
ered behind the integral sign for all possible powers satisfying YliiPi^Qi) = 
span an everywhere dense subspace, Tw^\ of the straightly defined com- 
plete function space, St("), which can be introduced by extending basis B 
into a basis B = {Qi, . . . , Qk} of the whole X-space and replacing the above 
complex polynomials by the real polynomials Y\{Qi: X)°'^ , where = n. 

(n) 

While depends on B, this straightly defined function space is in- 

dependent from the choice of the basis B. This function space naturally 




A{\X\,K)U^^\llzf^{Ku,X)zf{Ku,X))e'^'<'^>dK. (60) 
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emerges on the Cartesian product X x Z as Cartesian product of the corre- 
sponding function spaces. The same statements are true for the Z-Fourier 
transforms, ¥z{Tw^^) and ¥z{St^^^), of these function spaces. When these 
function spaces are defined for particular powers pi and qi, they are denoted 
by Fz(Tu;^* resp. Fz(§t(P-«»)). 

It is also observed in [Sz8] (cf. Theorem 6.1) that, for a function repre- 
sented by twisted Z-Fourier transform, j4(|X|,Er) is a uniquely determined 

L^-function. Thus, by fixing a basis B and two H-type groups H^""'^^ and 
H^'^ ^ satisfying {a + b) = {a' + 6'), there exists a well defined one to one 
correspondence kb '■ Fz(Twg^) ¥z{T'w^^) that maps a function ex- 
pressed by and complex structures J^^^ to functions where just 

the complex structures are exchanged for J^j'' ^ Let it be pointed out again 
that this map operates on functions defined by the twisted Z-Fourier trans- 
form and not on the ones standing behind the integral sign. This map has a 
unique extension defining a bijection between the straightly defined function 
spaces. 

In order to check out if this map is an intertwining operator, let the 
Laplacian be act on the twisted Z-Fourier transform formula. Like for Z- 
crystals, operator M = "^daD^* appears behind the integral sign in the 
form i\K\DK^*, which acts only on the polynomial part, resulting |Erjm, 
where m = ^{—pi + Qi)- Term involving Az appears inside as 47r^|-ftrp(l + 
Finally, the radial Laplacian in Ax acts on X-radial functions by 
radial differentiations and multiplications with radial functions, furthermore, 
the action of the X-spherical Laplacian is nothing but multiplication with 
the corresponding eigenvalue. Thus the above operator really intertwines 
the Laplacians (cf. these details in Section 8 of |Sz8| ). 

Although the considered functions are everywhere dense in the straightly 
defined complete function space, the boundary conditions can not be com- 
puted by them. The main problem is that, regarding the K-variable, func- 
tion A{\X\,K)uP([lzf'{Ku,X)-zl'{Ku,X)) is not the multiple of a sin- 
gle s*^-order spherical harmonic by a K-radial function even if one uses 
functions = (j){\X\,\K\)ip{Ku), where (p{Ky) is a spherical har- 

monic of degree / = deg{ip). This claim follows from the fact that func- 
tion Y\z^^{Ku-,X)'zf'{Ku-,X) is not derived from a homogeneous polyno- 
mial but, after performing the indicated powering imposed on functions 
Zi{Ku,X) = {Qi + Ufc^iQi), X) and its conjugate, it appears as a proper 
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sum of the form: 

llzr{K^,X)z^^'{K^,X)= (61) 

n 

(ai...a„) i a=J2ai=0 

where nj = {pi + Qi), n = X^rij, and index a = ^ Oj indicates that how 
many complex structures Jk^, appear in the term determined by the expo- 
nents (oi . . . a„). Sub-sum R^"'\Ku, X) consists exactly those terms where 
this number is a. Although these functions are derived from a-homogeneous 
functions, neither they nor ip{Ku)R^°'\Ku, X) are spherical harmonics re- 
garding Ku- One obtains the desired s^'^-order spherical harmonics by pro- 
jections and where compound index s = (s,/, a) indicates the 
degrees both of the target and original functions. Like the projections in- 
troduced regarding the X-variable, also the latter projections are the poly- 
nomials of the Laplacian A^^^. More precisely, they appear in the form 
^sII^^A^''^'^ where the first term is a constant and the last term pro- 
duces from a (/ + a)*'^-order polynomial of the K- variable an s^'^-order poly- 
nomial which is then projected to the s*'^-order spherical harmonics in the 
same way how it was defined regarding the X-variable. (These computations 
are described in Section 6.4 of |Sz8j .) 

The corresponding formula for constructing functions by which the bound- 
ary conditions can be handled is: 

/ Y,M\x\, lK|)n^(^(K.)n5)(n<'^f )(i^n,x))e'<^'^)dK. (62) 

It is also pointed out in |Sz8| that for fixed spherical harmonics ip{Ku), 
non-trivial projections are defined just for those s-values which satisfy the 
inequality (/ — a) < s < (/ + a) and s has the same parity as (/ — a) resp. 
(/ + a). Furthermore, functions defined by different a's are projected into 
independent subspaces. That is, this formula generically involves (a + 1)- 
tuples, (</>(/_a), . . . ,4>(fj^a))-> of functions depending \X\ and \K\. 

The functions constructed in this way form a larger space, F^(LTt(;g*''^'''), 
then those constructed in ()60p . Since functions ()60p form an everywhere 
dense subspace in the whole straightly defined space, they are everywhere 
dense also in space spanned by functions defined in (j62p . This means that in 
the newly defined space the functions can just be approximated by the pre- 
vious functions. However, no natural (obvious) approximation exist what 
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would make the consideration of the new formula evident. Anyhow, the 
statement about the density of these function spaces is enough to establish 
the following statement: For any two members of an isospectrality fam- 
ily, the intertwining operator defined for the primarily functions extends to 
the newly defined ones such that it associates functions defined in terms 
of same {X, Z) -radial functions and the corresponding projected functions. 
More precisely, this is the only option for a continuous extension. Indeed, if 
the operator is defined in this way, then on F^(Tti;g*''''^) it is the same as the 
operator defined originally. This argument implies that the operator defined 
in terms of the newly defined functions is still an operator intertwining the 
Laplacians. 

For investigating the boundary conditions there are two important tools 
established. One of them is the Hankel transform, which is proved in Section 
6.3 of |Sz8] in the following form: Any s^^ -order spherical harmonic Ci'^\Ku) 
defined on the unit sphere of defines, by the formula 

■Hl'\<P)i\Z\)ci'\z^) = [ cP{\K\)4'\K^)S^'^)dK, (63) 

a uniquely determined transformation 'H!f\(f)){\Z\) on the L?' -radial func- 
tions 0(1^1) which depends just on s and I. This statement implies that the 
newly defined function space is appropriate for constructing the complete 
space both of the Dirichlet and Z-Neumann functions in terms of the (X, Z)- 
radial functions. It turns out that these two conditions can be character- 
ized as being scalar, meaning that they are satisfied if and only if functions 
'h!'i^\(I)s){\X\, \Z\) appearing in the above formulas satisfy them individually. 
Actually, these functions are also explicitly described in [Sz8], establishing 
this statement a much stronger form. As a result, these functions are inter- 
twined by the above operator, indeed. 

The other important tool is the inner algorithm (cf. Section 6.7 in [Sz8]) 
induced by the action of the angular momentum operator or operator Dz* 
on the second Z-Fourier transform formula. In this algorithm this action is 
iterated such that in the A;*''-step the desired functions are constructed by 
functions (j)s obtained in the previous step by using Hankel transform, 
radial differentiation, and certain combinations of them which can be de- 
scribed as averaging by the roulette operator [SzSj . The point is that this 
process involves just the radial functions and it ends up either in finite or 
infinite steps. In the latter case the sought functions are obtained by the 
limiting A; — )• oo. 

There is far more difficult problem to establish the above statement also 
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for functions yielding the regular Neumann condition. The source of these 
difficulties is that, contrary to the Dirichlet and Z-Neumann conditions, this 
one does not break down to single individual functions, but, it can be ex- 
pressed in terms of all functions . . . More precisely, with 
the help of the inner algorithm, one can construct a compound Neumann 
operator Ni'^'"^((/>(j_a), . . . , acting on radial functions such that a 
function constructed by (j62p satisfies the regular Neumann condition if and 
only if the radial functions defined by the compound Neumann operator 
vanishes at the boundary points. That is, also this condition can be ex- 
pressed in terms of radial functions, therefore it is also intertwined by the 
above operator. 

The complexity of angular momentum operator is fascinating. Its 
action can be described by the Hankel transform and the above mentioned 
inner algorithm in a more precise way (cf. Sections 6.7 and 6.8 of [Sz8) ). 
These tools reveal that it appears as the sum of an extrinsic orbiting oper- 
ator, L, operating on the radial functions (ps without defining permutations 
(i. e., averaging by the roulette operator) on them, and an intrinsic spin 
operator, S, acting on radial functions, but contrary to L, it defines also 
an irreducible permutation on functions 0s . The actual appearance of this 
operator on the above functions is: 

/ Os(%-a),...,0{/+a))nJ^(^nJ)(n^f^O(i^«,^))e'<^'^>dK (64) 

The intrinsic life of particles is encoded into S. Also the strong nuclear 
forces, keeping the particles having interior together, can be explained by 
this operator. 

The rest part, (E = Ax + (1 + jl-'^P)^^ -|- L, of the complete Laplacian 
A exhibits just orbiting spin. Its action can be reduced to a single radial 
function. More precisely, the exterior operator CE, on constant radius Z-ball 
bundles reduces to a radial operator appearing in terms of the Dirichlet-, 
Neumann-, resp. mixed-condition-eigenvalues A,-''^ of the Z-ball Bz{R) in 
the form 

= Ati"{t) + {2k + An)i'{i) - (2m/i + + (65) 

By the substitution ^ = y xf''^ /4, this is exactly the radial Ginsburg- 
Landau-Zeeman operator (j58p obtained on Z-crystal models. 

The physical forces corresponding to operator CE are the weak nuclear 
forces by which the beta decays are explained. The theoretical establish- 
ment of this force went through enormous developments [Vl IWll IW2] . It 
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started out with Fermi's theory which was highly surpassed by the Glashow- 
Weinberg-Salam theory, whose greatest achievement was the unification of 
the weak nuclear force with the electromagnetic force. The above statement, 
which is unknown both in physics and mathematics, is an exact mathemat- 
ical establishment of this unification. In physics the unification with the 
other forces, that is with strong nuclear forces and gravitation, is one of the 
most intensely investigated unsolved questions. Since the Z-crystal and Z- 
ball-bundle models unify the electromagnetic and weak nuclear forces also 
with the strong nuclear forces, this unification is actually much stronger 
than those established in the Glashow-Weinberg-Salam theory. The main 
unifying idea is that all these forces can be derived from the very same oper- 
ator A. They are distinguished only by the invariant subspaces to which the 
A is restricted to. That is, the individual forces emerge on these individual 
invariant subspaces separately. 

Conclusion: On these mathematical-physical models the investigated isospec- 
tralities are equivalent to the spectral C/I-symmetry known in elementary 
particle physics. If the elements of the basis B are picked up from the in- 
variant subspaces X^"''^ and then the isospectrality is the manifestation 
of the pure spectral C-symmetry. In other words, the isospectrality proofs 
in my constructions mathematically demonstrate that, regarding the spec- 
trum of the Hamilton operator, these physical models obey the physical C- 
symmetry law. That is, this spectrum is not changing if some of the particles 
are exchanged for their antiparticles. 

7.4 Wave mechanics. 

The wave operators corresponding to the above Hamilton operators emerge 
in the Laplacian of the static resp. solvable extensions of the nilpotent 
groups. In physics, where this operator should be a hyperbolic wave opera- 
tor, the extended metric must be indefinite having Lorenz signature, where 
the time axis has signature —1. Regarding both bundles, the solutions of the 
wave equations considered on the whole bundle can be written up by explicit 
integral formulas (cf. Section 7 of |Sz8| ). On ball-bundles, corresponding to 
the static and solvable extensions, these de Broglie wave packets appear in 
the form 

[ 0,n(^Vn?^(n^r4O)e'^^'^'^^""*^rf^, resp. (66) 
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where w is a constant depending on the particle system. On Z-crystals, 
an appropriate discrete version of the integral defines these wave packets. 
Then, on the extended Z-crystals, the Laplacian is the sum of a Schrodinger 
operator determined for electron-positron systems and an electron-positron- 
neutrino operator accompanying the electron-positron system. On ball- 
bundles, the corresponding particles regarding the CE-operator are the W- 
and Z-particles, introduced by Weinberg in his weak-force theory. Analo- 
gous particles can be introduced regarding the complete Laplacian ffi + S. 
These details are completely described in Section 7 of |Sz8j . 

On manifolds endowed with indefinite metrics the isospectrality ques- 
tions can not be raised in general. However, this question can be raised re- 
garding positive definite extensions, where one extends the natural positive 
definite inner product defined at (0, 0, 1) on the tangent space. These man- 
ifolds can be considered as non-relativistic Newtonian space-time models. 
Actually the harmonic isospectral manifolds discussed in the AriasMarco- 
Schiith paper appear exactly among these manifolds. The metrics on nilpo- 
tent groups are neither harmonic nor Einstein thus they could not have been 
effected by the AM&S-theorem even if it was a right statement. The proof of 
isospectrality on the positive definite solvable extensions can be established 
similarly as on the nilpotent groups. In this case the radial functions should 
be considered in the form where t > is the new parame- 

ter added to (X, Z) which indicates that the solvable extensions are defined 
by the half-space extension of two-step nilpotent Lie groups. These details 
and the isospectralities on the boundaries of Z-ball-bundles are described in 
[Sil]. 

The problems arising in my papers [Sz4|, ISz5| . The intertwining 
operators are not correctly established in my papers cited above. These 
difficulties concern only these papers, where the isospectrality examples are 
constructed on Z-ball resp. Z-sphere bundles. They do not effect my other 
constructions performed on Z-torus bundles, and, even in these papers, they 
concern only the construction of the intertwining operators, while the other 
statements are unaffected. 

These operators are constructed in [Sz4llSz5] by means of a fixed complex 
structure Jq G Eskewi'^)- The complex structure obtained by a deformation 
is denoted by Jq. Then, the intertwining operator, k, is defined in these 
papers by the correspondence 
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k: A{\X\,Z)U^^\ll{Qi + iJo{Qi),Xy^{Q,-iJo{Qi),X)'^^) (67) 

i 

^ A{\XlZ)uP([l{Qi + i4{Qi),xr{Q,-i4{Qi),X)'^^) 

i 

where Qi ^ X are arbitrary X- vectors. That is, the associated functions 
appear in terms of the associated complex structures Jq resp. J'q in the very 
same form. It obviously intertwines these two complex structures 

It was H. Fiirstenau [F] who recognized that this map was not well 
defined. For if one considers a fixed complex basis B = {Bi, . . . , -6^/2} on the 
X-space, then the above correspondence can be prescribed only for functions 
defined by such Qj's which are in the real subspace SpanM.(B) spanned by 
real linear combinations of the basis-elements Bi. It is evident, that these 
correspondences will define the action of k on the rest of functions defined 
by other Qj's not lying in Span^CB). The very same map can be defined in 
terms of the complex valued polynomials constructed above regarding the 
basis B. That is, it is enough to consider only the correspondence which 
associate polynomials (appearing in the same form regarding the associated 
complex structures) to each other. It is also clear that the action of this 
well defined k on functions defined by Qi's not lying in Span^CB) differs 
from the action described in (I67p . Thus the above map is not well defined, 
indeed. (Fiirstenau used different arguments for explaining this problem.) 

This recognition triggered the reconstruction of the above ill-defined in- 
tertwining operator. Further investigations showed that a correct operator 
could have not been established by a single complex structure Jq. Instead, 
all complex structures Jz^ must be involved to its definition. Even the 
exchange of Jq for Jz^ does not alter this newly defined operator into a 
correct one. It becomes the desired intertwining operator, however, if one 
applies also the Z-Fourier transform to the latter altered functions. This 
step is really necessary because the operator defined in terms of functions 
standing behind the integral sign intertwines just those physical operators 
which appear after the Laplacian enters behind the integral sign. In terms 
of a complex basis B and polynomials written up regarding this basis, the 
reconstructed intertwining operator can be defined by the functions intro- 
duced by the twisted Z-Fourier transforms (j60p and (j62p . 

These reconstructed formulas were described first in |Sz6| and lectured 
about them also at the international conference held at CUNY, in February 
of 2006 |Sz7] . The physical content reviewed in this paper is the result of 
a recent development. It is described in |Sz8j . In this paper a completely 
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new mathematical idea, namely the Hankel transform, is introduced, by 
which both the electroweak and strong interaction Hamilton operators can 
be established. By this transform simple transparent proofs can be given also 
for mathematical theorems such as the density theorems and several other 
statements concerning intertwining of Laplacians or boundary conditions. 
In |Sz6| these statements are established by a different integral transform, 
called dual Radon transform. An other statement used there is the so called 
independence theorem. It should be pointed out that the exploration of 
the physical contents inherent in these structures is not quite much efficient 
by these tools than by the Hankel transform. This statement is certainly 
true for electroweak and strong interactions, which become understandable 
just by the Hankel transform. The strength of the latter method can be 
demonstrated also by the fact that all theorems established in [Sz6j can be 
established also by this tool. 

References 

[AM-S] T. Arias-Marco and D. Schiith: Local symmetry of harmonic spaces 
as determined by the spectra of small geodesic spheres. arXiv: 1001. 1611. 
Presented also at the International Conference on Spectral Geometry, 
July 19-23, 2010, Dartmouth College. 

[Be] A. L. Besse: Manifolds all of whose geodesies are closed. Ergebnisse 
der Math, und Grenz. 93, Springer Verlag, 1978. 

[BCG] G. Besson, G. Courtois, S. Gallot: Entropies et regidites des espaces 
localement symetriques de courbure strictement negative. Geom. Func. 
Anal.h: 731-799, 1995. 

[Bo] D. Bohm: Quantum Theory. Dover Publications, INC., New York, 
1979 

[DR] E. Damek and F. Ricci: A class of non-symmetric harmonic Rieman- 
nian spaces. Bull. Amer. Math. Soc. (NS. ):27; no. 1, 139-142, 1992. 

[F] H. Fiirstenau: Uber Isospektralitat von topologischen Ballen. Diploma 
Thesis Univ. Bonn, 2006. 

[GW] C. Gordon and E. Wilson: Continuous families of isospectral Rieman- 
nian metrics which are not locally isometric. J. Diff. Geom. 47: 504-529, 
1997. 



46 



[H] J. Heber: On harmonic and asymptotically harmonic homogeneous 
spaces. Geom. Fund. Anal.W: 869-890, 2006. 

[K] A. Kaplan: Riemannian nilmanifolds attached to Clifford modules. 
Geom. Dedicate, 11:127-136, 1981. 

[Kn] G. Knieper: New results on noncompact harmonic manifolds . 
\arXw:09T073F^ 

[LL] L. D. Landau, E. M. Lifshitz: Quantum Mechanics. Pergamon, 1958. 

[N] Y. Nikolayevsky: Two theorems on harmonic manifolds. Comm. Math. 
Helv.80: 29-50, 2005. 

[P] W. Pauli: Pauli Lectures on Physics. Dover, 2000. 

[RS] A. Ranjan and H. Shah: Harmonic manifolds with minimal horo- 
spheres. J. Geom. Anal.12: 683-694, 2002. 

[S] D. Schiith: Isospectral manifolds with different local geometries. J. 
Reine Angew. Math. 534: 41-94, 2001. 

[Szl] Z. I. Szabo: Lichnerowicz conjecture on harmonic manifolds, J. Diff. 
Geom. 31: 1-28, 1990. 

[Sz2] Z. I. Szabo: Spectral theory for operator families on Riemannian 
manifolds. Proc. of Symp. in Pure Math., 54: 615-665, 1993. 

[Sz3] Z. I. Szabo: Locally non-isometric yet super isospectral spaces. Geom. 
funct. anal. (GAFA), 9:185-214, 1999. 

[Sz4] Z. I. Szabo: Isospectral pairs of metrics on balls, spheres, and other 
manifolds with different local geometries. Ann. of Math., 154:437-475, 
2001. 

[Sz5] Z. I. Szabo: A cornucopia of isospectral pairs of metrics on spheres 
with different local geometries. Ann. of Math., 161:343-395, 2005. 

[Sz6] Z. L Szabo: Reconstructing the intertwining operator and new striking 
examples added to "Isospectral pairs of metrics on balls and spheres with 
different local geometries". DG/0510202. 

[Sz7] Z. I. Szabo: Laplace spectra on open and closed Zeeman manifolds. 
^aih:SP/0602474\ 



47 



[Sz8] Z. I. Szabo: Exact mathematical models of a unified quantum theory; 
Expanding and static micro universes, i i arXiv: 0909. 39(M 

[V] M. Veltman: Facts and Mysteries in Elementary Particle Physics. World 
Scientific, 2003. 

[Wl] S. Weinberg: Dreams of a Final Theory. Vintage Books, A Division 
of Random House, Inc., 1994. 

[W2] S. Weinberg: The Quantum Theory of Fields Cambridge Univ. 
Press, 1995. 



48 



